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O 1 La transformada de Laplace. 


Se entiende por transformada de Laplace £ una aplicación entre un cierto conjunto de 
funciones reales de variable real "x" definidas en [0, oo) y el conjunto de funciones de variable "p" 
compleja, p = o + i w, definidas en un semieje (a, co), o sea: para valores de o superiores a un cierto 
valor real a. Una característica principal de la transformada es la de ser una aplicación lineal entre 
espacios vectoriales. Posteriormente precisaremos ciertas condiciones suficientes para que las 
funciones originales f(x) tengan transformada F(p). El espacio imagen está caracterizado por funciones 
analíticas en el campo complejo pero nuestro interés excluye caracterizar tanto el dominio imagen de 
la transformada como el estudio de la imagen geométrica de la función racional F(p) en el plano de 
Gauss. En la mayor parte de los problemas y desde luego en los que estamos aquí interesados el que 
"p" sea considerada real o compleja no va a tener ninguna relevancia así que podrá ser visto como un 
simple parámetro —real, si se quiere— ya que la ventaja de considerar a "p" una variable compleja 
sobreviene casi exclusivamente al pretender disponer de una expresión general para la transformada 
inversa £”*. El tratamiento en el que estamos interesados implicará el manejo reiterativo de algunas 
tablas de transformadas bastante reducidas que nos va a permitir sistematizar racionalmente el uso de 
estos métodos de transformación aplicados al análisis de la flexión simple en barras de rigidez 
constante. 


$ 1.1 Algunas cuestiones referentes a la transformada £. 


La transformada de Laplace queda definida por la siguiente integral impropia: 
£ [£(9] = $, e P*1G9) dx = F(p) e 


que es, además una integral paramétrica dependiente de "p" y por lo tanto una función F(p) 
denominada "transformada de Laplace de f(x)".La aplicación £ está naturalmente restringida a las 
funciones que posean transformada lo que implica que las funciones transformables dependen de la 
convergencia de la integral “) a la que está supeditada la propia expresión de la transformada. 


El campo fundamental de las funciones que contemplamos preferentemente en los problemas 
técnicos unidimensionales es aquel que en un cierto intervalo [O, L] las funciones admiten en cada 
punto límites finitos a derecha e izquierda. Si estos límites en algún punto no son iguales habrá en ese 
punto una discontinuidad "de primera especie" —en el sentido clásico- o bien "un salto finito". 
Denominaremos a estas funciones "segmentariamente continuas" en el intervalo. Además, para 
garantizar la existencia de transformadas a estas funciones segmentariamente continuas se les deberá 
exigir la condición adicional de ser exponencialmente acotadas lo que significa que deben satisfacer 
para cada x > a el que existan dos números reales "C" y "a" tales que 


If) |<Ce**  obienl f(x) e**|<C 


El conjunto de las funciones segmentariamente continuas y exponencialmente acotadas 
forman un espacio vectorial real. Con f(x) perteneciente a este conjunto de funciones la integral 
US f(x) e PX dx es acotada y f(x) posee transformada de Laplace para p > a. Trivialmente si la 
función f(x) es acotada es también exponencialmente acotada para a. = O, así que si f(x) cumple el 


resto de condiciones admitirá transformada de Laplace para p > 0. Una cota de la integral Y será 
simplemente k/p siendo "k" la cota superior de f(x). 


Observamos, además, que el límite inferior de la integral *? es cero con lo que funciones como 
las representadas en la figura 1.1 siendo continuas en (0,00) tendrían la misma transformada. Para 
restablecer la uniformidad de la transformada inversa podemos considerar, en principio, que las 
funciones consideradas sean f(x) = 0 para x < 0, Posteriormente la función "escalón" —de Heaviside— 
nos permitirá resolver formalmente esta situación. 


2 2 
y=x+*a  0<x +a 0<x e 


y = xa x<0 
figura 1.1 


Resumiendo: podemos exigir como condiciones suficientes para que una función f(x) tenga 
transformada de Laplace las siguientes 


- que para x<0 sea f(x) = 0 

- que en [0,00) sea segmentariamente continua pudiendo existir un número 
finito de puntos de discontinuidad de primera especie. 

- que en [0,00) sea exponencialmente acotada. 


Estas precisiones matemáticas no van a interferir en nuestro campo de trabajo. La mayoría de 
las funciones utilizadas normalmente en la técnica y, desde luego, las que aquí vamos a considerar 
poseen transformada y esto es consecuencia, entre otras cosas, de que dichas funciones decrecen 
menos rápidamente que la función exponencial e ”* teniendo en cuenta que podemos elegir el valor 
"p" de dicha función. Esta elección del "p" suficiente, incide, como ya se ha sugerido, en el intervalo 
del plano imagen en el que F(p) queda bien definida. Esta última cuestión puede ponerse de manifiesto 
si el limx30 € P* f(x) = 0 para valores convenientes de p> a. 


Finalmente podemos añadir que desde el punto de vista operacional la transformada £ puede 
ser considerada un operador integral caracterizado desde la forma general 


E(p) =$, np, 0 £G9 dx 


donde n(p, x) es una función analítica que en el caso de la transformada de Laplace es la función e P* 
siendo el intervalo m = [0, oo). Otras transformadas notables son las de Fourier, las finitas de Fourier 
sen y cos, la de Mellin, la de Hankel o la de Legendre...etc. Junto a ellas puede recordarse la 
transformada 3 definida por 3 [£()] = Ni=of() z* que ocupa, en cierto modo, el mismo papel que £ 
para el estudio de sistemas con señales discretas. 


$8 1.2 Campo de aplicación de la transformada. 


El contexto más idóneo de la transformada de Laplace resulta ser el de las ecuaciones 
diferenciales consecuencia de modelos y dispositivos físicos de respuesta lineal. En los sistemas 
dinámicos la variable 'x" es la temporal con lo cual que se comprende que físicamente la restricción 
f(x) = 0 para x< 0 resulta obvia. Asimismo, en los problemas mecánicos o elásticos en los que estamos 
interesados como el de la viga-ménsula de sección constante y longitud "L" de la figura 1.2, una vez 
asimilada a un modelo unidimensional, la deformación y(x) en el intervalo [0, L], carece de sentido 
para x<0 así como para x > L. 


y(x) 


figura 1.2 


La utilización usual de la transformada se presenta en la resolución de problemas regidos por 
ecuaciones diferenciales lineales o integro-diferenciales en el dominio de la variable "x". Mediante la 
transformada £ llevamos esta ecuación al dominio de la variable "p" donde se resuelve con 


operaciones algebraicas elementales y posteriormente con la transformada inversa £”* se restituye al 
dominio original según se muestra en el diagrama del cuadro 1. 


condiciones inciales 


ecuación diferencial ordinaria 
o 


integro-diferencial en "x" ecuación algebraica en "p" 


en 


á 


solución general cálculos 


solución en "x" solución en "p" 


CUADRO 1 


El paso entre ambos dominios se realiza en los casos más sencillos mediante tablas con 
operaciones complementarias de operaciones para su aprovechamiento exhaustivo. En el caso de 
ecuaciones diferenciales ordinarias lineales con condiciones iniciales dadas, esto es: en los problemas 
"de valor inicial", la resolución presenta bastantes ventajas sobre los métodos clásicos hasta el punto 
de que el estudio de un amplio conjunto de sistemas lineales eléctricos, electromecánicos, de control 
etc. se precisan y se analizan en el propio dominio de la "p" (funciones de transferencia, mapas de 
polos y ceros, etc.). Estas ventajas se pueden poner de manifiesto, en una primera estimación, en los 
siguientes aspectos: 


- La solución puede ser obtenida de forma sistemática mediante tablas, en este sentido y 
salvando las distancias se asemejaría a la antigua utilización de logaritmos. 

- Las condiciones iniciales se incorporan a la propia transformación. No hay determinación 
posterior de constantes. 


- En determinados problemas, la transformada puede ofrecer información más amplia y 
completa en el propio dominio de la "p" 


En el análisis de la flexión simple de barras rectas una vez reducido el problema a un 
modelo unidimensional lineal con condiciones en los extremos, tendremos que adaptar los 
métodos naturales de la transformada para su utilización sistemática a problemas que no son 
problemas "de condiciones iniciales" sino de condiciones en los extremos o "de contorno”. 
Esta adaptación se realizará en el O 5 estableciendo previamente los distintos tipos estándar 
de sustentación en las vigas. 


Ahora, nuestro primer paso será construir esas tablas de transformadas tanto para las funciones 
elementales que centran nuestro interés como para los operadores. 


8 1.3 Transformadas de los operadores 
$ 1.3.1 Derivación e Integración. 


Suponemos conocida £ [f(x)] = F(p) y queremos hallar la transformada de la derivada siendo 
f(x) diferenciable en x > 0 y f (x) transformable Laplace. Aplicando la definición 


£ [69] = f, e P*f'Go dx 

y resolviendo esta integral por partes, obtenemos 

£ [F69] =P £60| A + Só p e PX f(x) dx 

como debe ser el lim, e P* f(x) = 0, resulta 

£ [£(9] = —£(0) + p $ e P*£G9 dx = p E(p) — £(0) 


Este resultado puede ser reiterado para hallar f(x) supuesta f(x) derivable dos veces en x > O y 
f"() transformable. Siendo g(x) = f'(x) 


£ [19] = p £ [£'G)] — £(0) = p (p F(p) — f(0)) — £'(0) =p F(p) — p f(0) — f1(0) 
y, supuesta f(x) diferenciable "n" veces en x > 0 y la existencia de £ [f "(x)], se obtendrá 


£ [E 69] = p* F(p) =p" * 10) =p" *£ (0) —...— PE %0) -£%%0) 


Utilizaremos para la función y(x) las siguientes transformaciones, con £[y(x)] = Y(p): 


y Oia b o D V(D) — y(0); y 0) E Ear e Y (p) —- p y(0) — py (0) — y"(0) 
y"(%).....£.....p" Y(p) — P y(0) — y'(0)  y"G0.....£.....p* Y(p) — p” y(0) — p” y '(0) — p y"(0) — y""(0) 
En cada caso hemos considerado la continuidad de las funciones f(x), f'Gx),..., f "() para 


valores x > 0O.Estamos, sin embargo, tratando con funciones segmentariamente continuas que sean 
transformables así que cuando aparezcan discontinuidades de salto finito estas fórmulas deberán ser 
corregidas. Estas correcciones se explicarán más adelante y se superarán formalmente con el auxilio 
de las funciones "escalón" e "impulso" por lo que estas precisiones pueden no ser tenidas ahora en 
cuenta. No obstante si f(x) es discontinua en x = a punto de encuentro de sectores continuos con 
límites finitos distintos a la derecha f(a”) y a la izquierda f(a”) del punto en cuestión —figura 1.3-, la 
fórmula correcta de aplicación de la transformada de la primera derivada sería 


£ [£'(x)] = p F(p) — £(0) - e **(f(a”) —f(a)) 


De la misma manera, si es f(x) diferenciable en x > 0 pero f'(x) no es continua en x = a —con 
salto finito—, siendo f"(x) transformable, obtendríamos la fórmula 


€ [£"(x9] =p E(p) — p £(0) —£1(0) - e *>(£ (a) -f' (a) etc... 


En cada caso se observa que el factor g(a*) —g(a”) es la contribución de las discontinuidades, 
siendo su valor absoluto el denominado "salto" de la función discontinua. En nuestro contexto estas 
discontinuidades estarán generalmente formadas por dos funciones diferentes a derecha e izquierda 


cuyos valores en el punto son distintos con lo que el salto será de la forma f;+: (ax) — f;-1(ax) tal y como 
se muestra en la figura 1.4. 


figura 1.3 


Supuesto ahora f(x) integrable en x >0 y dE f(s) ds una función transformable, la 
transformada de la integral: 


£ [f, £(s) ds] =S, e PX f¿F(S) ds) dx 


puede ser realizada por partes obteniendo 


£ [9 169 dx] == 5 e7P% ¿6 del y +77 er f09 dx 


Si f(x) tiene transformada limy-,0 e P* UN f(s) ds = 0 y de aquí se deduce que 
x _ 1d po _ F(p) 

£ [f, £G0 dx] e dE de A 

y se puede escribir además 


£ [f, dx S, dx.....(k veces)..... ES) dx] = AS 


y 


figura 1.4 


$ 1.3.2 Traslación en el dominio de la "x" y en el de la "p". 


Conocida £ [f(x)] = F(p), la transformada de £ [f(x—a)] será la integral 
£ [f (x—a)] = fee f(x-a) dx=e P? S. eP (3) f(x-— a) d(x-a) = e P? E(p) 


Este resultado es válido solamente para x > a. Sia la función f(x) le hemos exigido que sea f(x) = 0 
para x< O deberá ser f(x — a) = O para x< a, como se muestra en la figura 1.5, de manera que el 
resultado sería válido en todo el semieje real positivo. Podemos presentar, no obstante, la transformada 
en la forma 


£ [f (x —a)] = e P? F(p), a>0, siendo f(x — a) = 0 parax<a 


y f(x - a) parax>a 


A 


y 


f(x) para x> 0 
[Oparax<0] 


figura 1.5 


Como en otras transformadas, este resultado será formalmente enunciado sin ninguna 
precisión adicional, de manera muy simple con auxilio de la función escalón en el Dd 2. 
Análogamente podremos escribir 


£ [ft9 e **] = 1 e PX AX f(x) dx = Ie e” P+)x f(x) dx = E(p+a) con a> 0 


Ambos resultados: traslación en el dominio de la "x" y en el de la "p", que nos permitirán 
generar nuevas transformadas, se garantizan como máximo en el semiplano complejo derecho de la 
"p" o a partir de p > 0 si f(x) es acotada. 


$ 1.3.3 Transformada de un producto de funciones. 


Hallar £ [£1) £(9)] disponiendo de £ [f:(x)] = Fi(p) y £ [£269] = F2(p) requiere del 
conocimiento de la integración en el campo complejo y es consecuencia de abordar el estudio general 
de la transformada inversa £”*. Su aplicación estricta no es frecuente en el tipo de problemas en el que 
estamos interesados por lo que quedaría aquí fuera de contexto. No obstante, es posible encontrar 
dicha transformada por operaciones elementales. Así ocurre cuando la transformada racional de al 
menos una de las funciones, por ejemplo f,(x) pueda ser descompuesta en fracciones y no posea más 
que polos simples; es decir si 

B C 


E 
AL ei EE 


Sabiendo que la transformada de una constante es £ [k] = . (resultado que se verá posteriormente en 
el $ 1.4.1), teniendo en cuenta que £ [f(x) e **] = F(p + a) y la linealidad de £, entonces f,(x) será: 
(69 =A E *M+BeP+Ce +, 

y, en este caso, se obtiene de forma inmediata: 

£ [60 £,60] = £[Ae ** 60) + Be?* f(x) + Ce * Lx) +...]J=A Fo(p-a) + B Fa(p-b) + C Fa(p—c) +... 


sea cual sea f(x) a condición de que posea transformada F»(p). Esta transformación es representativa 
de algunas operaciones usuales en el dominio de la "p". Posteriormente podremos ampliar esta 
cuestión con un carácter más general. 


$8 1.3.4 Producto de transformadas. La integral de la convolución. 


Uno de los resultados más importantes de la teoría de transformadas de Laplace aparece al 
estudiar la función primitiva de un producto de dos funciones transformadas F;(p) y F2(p) en función 
de sus originales. El resultado obtenido no jugará solamente un papel en el cálculo sino que puede 
configurar todo un desarrollo en cierto modo independiente y dando lugar a una nueva operación que 
se proyecta en muchos problemas físicos. Se trata de hallar £”* [F,(p) F2(p)] siendo £[f:(x)] = Fi(p) y 
£[£09] = Exp). 


Podemos expresar dicho producto en la forma: 
Ex(p) Fa(p) = $, e7P" fu) du $ e-PYf(v) dv 


y transformar este producto en una integral doble en el recinto R(u, v): 0< u< oo, O<v< «o —según la 
figura 1.6 (a— con lo que este producto quedaría 


Esp) Ftp) = Suso Syzo PY (a) £a(v) du dv= SJ, y 1 (0) E (v) e P+D du dv 


Este desdoblamiento de las variables antecede para nuestros propósitos a un plegamiento del recinto 
mediante el cambio de variable definido por u = x— T , V= 1, aplicación biyectiva y bicontinua, 
continuamente diferenciable, siendo el Jacobiano de la transformación igual a la unidad, así que 
podemos pasar de las coordenadas u,v a las x,t y el cambio da lugar a la integral: 


Es(p) Exp) = f ¡e ePXf, (x—0) £, (1) dx dr 


el recinto u > 0, v > 0 se transformará en x > 1, T > 0 que es el representado en la figura 1.6 (b 


figura 1.6 


Podemos expresar ahora el producto de transformadas en la forma: 
Esp) FAP) = So Sap e PEE 0) ata) dx de = Je PX ( SE (=D fo(1) de) ax 
así que el resultado buscado puede venir expresado mediante una integral 
EL FP) Fap)] =$, £: 00) £a(1) de 


que es una integral dependiente del parámetro "x" y por ello la función de "x" buscada. Esta igualdad 
puede ser vista en la forma directa: 


EL. £, G:—1) £2 (1) d1] = F(p) FAp) 


Puede ahora formalizarse una operación denominada "producto de convolución" o 
simplemente "convolución"” de las funciones f(x) y f,(x) 


1169 + 69) = f, 60 £ (1) dr 
y expresar el resultado de la transformación £Tf,(x) * £269] = F(p) Fx(p) = £[£,609] £[£2G01. 


La operación así definida puede ser estudiada para averiguar la estructura que aporta a cierta 
clase de funciones. Si las dos funciones son continuas su producto-convolución existe. En estas 
condiciones esta operación resulta cerrada, conmutativa, asociativa y distributiva. Para que no siendo 


continua las funciones originales la integral sea convergente será suficiente que ambas sean acotadas 
en Cualquier intervalo del semieje real positivo. Proporcionar una interpretación "física" de la 
convolución nos llevaría a plantear, fuera de nuestro contexto, problemas muy particularizados. 


$ 1.3.5 Cambio de escala en "x" y en "p". 

Conocida la transformada F(p) de la función f(x), la transformada de fía x) y f£( A ) serán 
£ [f(ax)] = fe Pr f(ax) dx == Feos “Xfax) d(ax) == F(É) , az0 
El) = Ser ax=a ff Pa 1) ad) =aF (ap) , az0 

8 1.3.6 Multiplicación y división por "x". 

La transformada de x f(x) puede resolverse derivado la igualdad 

£ [109] = fe” P*£G0 dx =F(p) 

y se obtendrá 

4FCp) = 549 oe EG) dx= Pe PXxfG0 dx =- £ [xfG9)] 


dp 0 


_ a F(p) 


así que £ [x f(x)] = ab 


Teniendo en cuenta las sucesivas derivaciones respecto del parámetro "p" de los productos 
e PX f(x), podemos generalizar el resultado anterior para cualquier x” f(x) con "n" entero positivo, 
obteniendo sin dificultad 


e 691 =C1" A Y" E) 

De forma análoga, el cálculo de la transformada de 2 conocida la de f(x) se resolverá integrando 
adecuadamente en "p" la igualdad 

Se P% £(x) dx = E(p) 

es decir, escribiendo 

ES dk f, e *f(9 dx = $ FO) dk 

que puede realizarse en la siguiente forma 


E SE dra 10) Ez cala ea 


así que £ [2] = le E(p) dp 


ES 
$ 1.4 Transformadas de las funciones. 


Se trata ahora de encontrar las transformadas de un pequeño número de funciones elementales 
que son las que normalmente aparecen en las aplicaciones técnicas: polinómicas, exponenciales, 
circulares e hiperbólicas etc. que junto a la tabla de transformadas operacionales vistas anteriormente: 
cambios de escala en "x” y en "p", multiplicación por x” con "n" entero, producto por la exponencial 
etc. nos permitan generar un más amplio conjunto de transformadas con menor cantidad de 


información inicial. 
$ 1.4.1 Transformada de una constante. 


Siendo k un número real, la transformada 


a MP 
£1=J, Pdo kero 


En el campo de las transformadas de Laplace este resultado es algo impreciso: "k" debería ser un 
ai 27 A 1 B 
escalar que multiplique a una función nula para x < O que tendría por transformada E Una vez más, 


esta función será la función "escalón" unitario u(x) que estudiaremos en el $ 2.1 junto a algunas 
precisiones sobre su uso. Veremos que esta función escalón unitario vendrá representada por u(x) y 


que £ [u(x)] = a por lo que con mayor precisión escribiremos: 
£ [ku()] =k£ [ubo] = . 


esto garantizará, entre otras cosas, que la función original es nula para x < O como veremos 
posteriormente. En cualquier caso, el papel jugado por u(x) pasará a ser parte muy activa en el texto. 
Dicho lo cual, podremos prescindir de esta función en este y en otros casos prácticos cuando no sea 
imprescindible o cuando solamente adquiera un papel puramente formal en el cálculo como se sugirió 
en el $ 1.2. 


$8 1.4.2 Transformadas de las funciones potenciales. Transformada de un polinomio. 
Siendo "a" un número real positivo, la transformada de x* será 
a-(“.pxya 
£L]=/ o EP xtdx 
haciendo el cambio de variable p x = y, quedará para a > 0 


a 0 e cd 0% 1 E 1 
£lx1=5 (Y) e e CEN 


siendo I” la función "gamma" denominada también integral euleriana de 2* especie, convergente para 
valores mayores que cero y cuya primitiva no es expresable para valores no enteros mediante otras 
"7" 


funciones que no sean la propia Il”. Cuando "a" es entero ['(a) es el factorial (a —1)!. Cuando "a" es no 
entero se descompone usualmente en la parte entera más el resto de manera que solamente queda por 


utilizar una tabla de I' para valores entre 0 y 1. En particular siendo, por ejemplo, I( 5) = Vr, será 


Nuestro interés se centra en exponentes "n" enteros positivos para los cuales 


n! 
pn+1 


£ [x”] = con n entero > 0 


Serán de gran utilidad en las aplicaciones las transformadas particularizadas: 


Jl 
xX E REO IEC p? 
2 
x A e e e DU, p? 
6 
Xx aii p* 


En consecuencia, teniendo en cuenta lo dicho, la transformada de un polinomio con 
coeficientes reales P,(x) = a, x + ay xXx 4+ ....+ 21 X + ap, será 
_ apr! ar: (r-1)! as 
E POIS o A ta o 
Este resultado nos permite encontrar las transformadas de desarrollos en serie entera, 


e, LEO) a 
f(x) = 2n=0 X 


es s =x?2 
que nos pueden ser útiles para obtener las transformadas de funciones tales como a%,e”*%”...etc. 


Bastará tener en cuenta la transformada 


(0) m_f*(0) 
n! Al pn+1 


£l 


y aplicar el carácter lineal de £. 
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8 1.4.3 Transformadas de las funciones exponenciales. 


Aplicando los resultados obtenidos anteriormente para la transformada de una constante y la traslación 
en "p" 


e**f(x) A 


Ha] 

= 
ES 
[ab] 

Y 


podemos obtener directamente la transformada de la función e ** 
1 


£ [e **]==— con a>0 
p+a 


y para la transformada de e ** podemos proceder de la misma forma 
ca ERE 
£ [e*? ] = 24 


resultado válido solamente para p>a>0 


$8 1.4.4 Transformadas de las funciones trigonométricas sen Wx y COS WX. 


Para obtener las transformadas de las funciones trigonométricas sen wx y Cos wx puede 
procederse sobre la definición con una doble integración por partes o bien utilizar 


1 
e”* O E RR p-a 
aplicándolo a la función e'**= cos wx + i sen wx, es decir utilizando 
1) _ p+iw 
mx £ p-iw  p2+w2 
igualando las partes real e imaginaria se obtienen directamente 
w 
£ [cos wx ] = ; £ [sen wx ] = =— 
[ ] p2+w2 a [ ] p?+w2 


$ 1.5 Generación de otras transformadas. 


Las transformadas de las operaciones más significativas junto las transformadas de algunas de 
las funciones elementales nos permiten obtener nuevas transformadas sin tener que resolver las 
integraciones implicadas en la definición. Para la mayoría de funciones que utilizamos habitualmente 
en la técnica bastarán dos tablas muy económicas; una de funciones y otra de operaciones. Así, por 
ejemplo, teniendo en cuenta 

dF a 1 
Aids ista O A A econ ct do E 
dp p+a 
la transformada de x e ** deberá ser: 
1 
da 
e + 1 

£ [xe**]=-2= 
dp (p+a)? 


resultado que podemos obtener también utilizando las trasformaciones 


¡IMA A a AAA il 
y con alcance general para "n" entero positivo 
n ,-ax73-_ n! 
£lx'e ds (p+a)n+1 


Teniendo en cuenta la transformada f(x) e ** y las transformadas de sen wx y Cos wx 
halladas anteriormente, obtenemos también directamente: 
w 
(p+a)2+w2 


p+a 


' ax = 
: £ [e ** cos wx] PATO 


£ [e ** sen wx] = 
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y considerando, además 


d F(p) 
O Eoiiantón. dp 
resultan las siguientes transformadas 
-_2Pp y - _p?w 
£ [x sen wx] = e24w2y? : £ [x cos wx] = o21w2y? 


Finalmente teniendo en cuenta las transformadas anteriores se pueden obtener fácilmente: 


2 (p+a) 


l (p+a)?w? 
(p+2)24+w2)? ? 


£ [xe ** cos wx] = rnzrw2y 


£ [xe ** sen wx] = 


$ 1.6 Teoremas del "valor inicial” y del "valor final". 


Si f(x) tiene transformada F(p) y es f(x) transformable, podemos obtener dos resultados 
importantes que pueden deducirse de la expresión de la transformada de la derivada 


£ [£(9]= f, e P*f'G9 dx = p E(p) — £(0) 
y puesto que f'(x) tiene transformada, el lim, S. e PXf'(x) dx=0, así que 


0 lim: p E(p) 0 £(0) 
y, de aquí, supuesta la continuidad de f(x) en x = 0 
limy,0 f(x) = limp 320 Pp F(p) 


resultado conocido como el "teorema del valor inicial" 


Si en la expresión de la transformada de la derivada calculamos el límite para p > 0 

limp > 0 $, e PXf'G9 dx = f,"f'Go) dx = f(9]| 07 lim f(x) —f(0) = limp_,0 p F(p) —f(0) 
X>0%0 

habiendo supuesto la continuidad de f(x) en x = 0 de valor f(0) y por tanto 
liMx 700 f(x) 5 limo p E(p) 
denominado "teorema del valor final" dónde han de asegurarse la existencia de los límites implicados 
en la igualdad. Una forma de garantizar la última fórmula es que los polos del producto p F(p) se 
encuentren en la mitad izquierda del plano complejo de "p" ya que la presencia de polos en la mitad 
derecha corresponde a una función f(x) que crece —exponencialmente— sin límite cuando x crece y la 
presencia de polos sobre el eje imaginario significa la contribución de una oscilación permanente en 


f(x). Se trata de garantizar que p F(p) sea analítica en la mitad derecha del plano complejo de "p" 
incluyendo el eje o =0 de p=0+iw. 


$ 1.7 La transformada inversa. 


Aquí se considera esencialmente la aplicación de la transformada a los problemas regidos por 
ecuaciones diferenciales lineales con coeficientes constantes y en estos casos F(p) resulta una función 
racional en general irreductible. Una utilización sistemática de las tablas de transformadas de las 
funciones y de los operadores es suficiente para encontrar las funciones originales en el dominio de la 
"x" y, en general, una tabla reducida a lo indispensable de transformadas inversas puede facilitar la 
lectura inmediata de las mismas para este tipo de funciones. En general, la obtención de transformadas 
inversas puede abordarse con diferentes alternativas que son complementarias y compatibles entre sí. 
Aquí bastará describir este cálculo mediante la descomposición previa de F(p) racional en sumas de 
fracciones simples o, en su caso, por desarrollos de F(p) en serie y aprovechar los resultados de las 
transformadas directas obtenidas para las funciones y para las operaciones vistas anteriormente 
incluyendo la del producto de transformadas. 


A estos procedimientos habría que añadir algunas notas sobre los mapas de polos y ceros de la 
función F(p) que tienen cierto interés en determinado tipo de problemas físicos. Finalmente puede 
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darse alguna referencia al llamado método "de los residuos” semejante desde un cierto punto de vista 
práctico al de la descomposición en fracciones simples. 


$ 1.7.1 Desarrollo en fracciones simples. 


Siendo £”* una transformación asimismo lineal, si abordamos la transformada inversa de la 
fracción irreducible F(p) que no esté directamente en una tabla reducida de transformadas el paso más 
inmediato es el de convertir esta función racional en suma de fracciones simples que sí encontraremos 
en dicha tabla. 


Consideramos la función F(p) = EP) donde R(p) y Q(p) son polinomios de coeficientes reales 


Q(p) 
siendo el grado de R(p) menor que el de Q(p). Necesitaremos el conocimiento de los ceros de Q(p) 


que son las raíces de Q(p) = 0 y que serán los polos "p;" de F(p). Supuestos estos valores reales, la 
función racional resultará ser una suma de la contribución de cada p; con orden de multiplicidad "k" 

PEPA 

ppp (py? (p-p)?  P-Pi 


Pisa 


así que bastará con disponer en una tabla de transformadas inversas de funciones Gras 


y para estas funciones es suficiente tener en cuenta 
n! 
A pato; (EA. es Pp) 


de las que se deduce fácilmente 


A n1 ,-ax3- A 
ls 4 IS as 


. E ps A ; 
Cuando el polo "p;," es simple el numerador correspondiente a la fracción 5-p, “on la que contribuye a 


1 
la descomposición en fracciones simples puede venir expresado por la siguiente fórmula 
— R(pi) - R(p) 
Qi) QQ p= Pi 
y esta es la alternativa que permite reflejar la "primera regla de Heaviside" o "teorema del desarrollo" 
que se enuncia de la siguiente manera: 


"Siendo F(p) TO) , siendo el grado de R(p) es menor que el de Q(p), si todas las raíces de Q(p) = 0 


son reales distintas Py, PD», P3,...,Pn, entonces £” [F(p)] = 


p; raíz simple de Q(p)=0 


n R(D; ) PjX " 
-, —— ePj 
J=1 0 1(p;) 


Cuando existan pares de raíces complejas conjugadas a + i b, su contribución individual puede 
ser representada por la secuencia anterior con coeficientes complejos pero, como se sabe, se prefiere 
directamente describir la contribución de cada par a +1 b en el campo real por la fracción 

Mp+N 
—— , conM y N reales. 
(p-a)2+b2 


Xx 


Partiendo de las transformadas de las funciones e€e** sen wx ,e** cos wx vistas 
anteriormente y teniendo en cuenta el carácter lineal de £ se llega con facilidad a 


e a SU 
£ [Me ** (cos bx — £ sen bx)] = TEE 
y de aquí a 
ar N-Ma — _MPp+N 
£ [e ** (M cos bx + sendx)l= yo? 


En determinados problemas se prefiere escribir para las funciones en "x" anteriores la forma 
K e ** sen (bx + q) con K = : yJ(N -M a)? + M?b? , tg4= —. 


Análogamente puede procederse con las raíces complejas múltiples. 
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Anotamos aquí como un recurso adicional el de los desarrollos en serie de las funciones F(p) y 
la posterior transformación inversa de los términos de este desarrollo. Para estos desarrollos es 
fundamentar el contar con la transformación 

Anal A 


Básicamente esta transformada es la que justifica la "segunda regla de Heaviside" que puede 
transcribirse de la siguiente manera: 


"Siendo la serie F(p) = X5= Me pa Puna serie de potencias convergente de p, se verifica que 


-1 = An n-1n 
£ [E(p)] => n= 0 (m- 1)! Xx 


$8 1.7.2 Los mapas de polos-ceros. 


Sin entrar en consideraciones que nos alejan de la descripción práctica de la transformada y 
conservando a "p" simplemente como un parámetro cuyos valores pueden ser representados en el 
R(p) 


campo complejo, el caso que nos ocupa es el de las fracciones irreducibles de la forma E(p) = OY 


Recordando que los polos son simplemente las raíces "p" de la ecuación Q(p) = 0 y visto todo lo 
anterior, se comprende que la forma de la función original f(x) = £”'[ F(p)] queda determinada por 
estas raíces. El numerador no tiene influencia sobre la forma de la función original y por lo tanto no la 
tienen los ceros de F(p) que son las raíces de R(p) = O. Estos últimos afectarían solamente a la 
amplitud y al desfase de una contribución del tipo A e”** sen (bx + q). Esta dependencia de f(x) de la 
situación de los polos de F(p) nos lleva a representar sobre el plano de Gauss su "mapa de polos- 
ceros". Con esta imagen geométrica no solamente puede reconstruirse la propia F(p) sino también 
deducir su transformada inversa por métodos gráficos. 


A 
(p+1)? (p+2) p 

sobre el plano complejo será los de la figura 1.7 (a, donde los polos se representan por una "x" en 
tanto "o" será los ceros y entre paréntesis se precisa el orden de multiplicidad de cada raíz. El 


Tomando por ejemplo la función F(p) = , la representación de sus polos y ceros 


% “e, _ p+3 - A 
diagrama de la función F(p) = Sade ía el representado en la figura 1.7 (b. 
wW w 
(2) 
o o 
a figura 1.7 (b 
En los sistema lineales dinámicos en los que la variable "x" es la temporal, el concepto de 


estabilidad es esencial y en estos casos la existencia de polos en la mitad derecha del plano de 
representación proporciona sistemas que crecen con el tiempo por lo que el sistema resulta inestable, 
los polos de la mitad negativa del eje real contribuyen con términos t”* e” ** según el orden de 
multiplicidad "n" del polo. Los polos complejos conjugados en la mitad izquierda del plano 
proporcionan Pe exponencialmente amortiguadas, los simples imaginarios puros proporcionan 
ondas senoidales sostenidas, los múltiples imaginarios puros proporcionan términos crecientes no 
exponenciales y los polos en la parte derecha del plano proporcionan términos crecientes en forma 
exponencial o exponencial-ondulatoria (figura 1.8). 
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(n) 


(n) 


(a (e 


figura 1.8 


8 1.7.3 Acerca de la fórmula de Mellin-Fourier y del método de los residuos. 


Bajo determinadas condiciones, es posible obtener de forma general la función f(x) 
transformada inversa de F(p).Estas condiciones están en relación directa con las exigidas a la función 
f(x) para que tenga imagen en el dominio de la "p". La fórmula que permite obtenerla es conocida 
como "fórmula de Mellin-Fourier” 


Ep = 4 


2niv01-10 


APR) ex dz =f(x) ,x>0 

con z = 0 +iw, variable compleja (o y w reales). Es una integral a lo largo de la recta o = o; en el 
plano de Gauss, integral que depende del parámetro real "x". El valor o, es 0, > a, donde a es el valor 
que aseguraba la convergencia y por tanto la existencia de la integral S, GO! e” PX dx. La condición 
esencial para que esta fórmula sea válida es que el itinerario al que se extiende la integral, esto es: la 
recta perpendicular al eje real por 6, quede a la derecha de todas las posibles singularidades de la 
función F(p) y esto significa que F(p) es analítica en el semiplano a la derecha de dicha recta. Elegido 


01 de esta manera £”* [F(p)] = 0 para x < O que es la condición impuesta en un principio para la que la 
función tenga transformada. 


contorno Y 
o01—iw 


figura 1.9 


Disponiendo de esta fórmula para la transformada inversa, se trata ahora de calcular de la 
forma más sencilla la integral del segundo miembro. Consideremos ahora la integral 


Sy F(p)eP* dp ,x>0 


extendida al camino cerrado "W" de la figura 1.9 conocido como "contorno de Bromwich" que puede 
ser expresada como suma de dos integrales 


1 
211 
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A Ep) eP* dp += f, F(p) eP* dp 


2rmi01-10 


donde la curva y es una semicircunferencia de radio R >00 a la izquierda de la recta o = 0. Si F(p) es 
una función racional irreducible, entonces para x > 0 


; 1 

limz-,0773 4, FG) eP* dp=0 

resultado que es una aplicación directa e inmediata del lema de Jordan. Así que para fracciones 
racionales irreducibles: 

EE) = SI Fa) ez* dz == f,, F(p) eP* dp 


01-10 


siempre que x > 0 y o = o; deje a su izquierda a todas las singularidades de F(p). Esto nos permite 
calcular la transformada inversa por medio de la integral extendida al contorno cerrado '"W" que, en 
las condiciones dadas, puede evaluarse por medio del cálculo de los residuos alrededor de los polos de 


de ePX ya que 
ET — fa F(p) eP* dp = Zextendida resíduos de F(p) eP* 


todos 
los polos 


y, en consecuencia £”* [E(p)] = A co Sk G paa = ((p- p;)” F(p) eP* | el 


La sumatoria está extendida a todos los polos p; de F(p) con órdenes de multiplicidad "n" y hay una 
contribución a la sumatoria por cada polo. 


Aunque desde ámbitos totalmente diferentes, se puede observar la semejanza entre este 
procedimiento y el desarrollo en fracciones simples visto en el 8 1.6.1. En general, ambos métodos son 
muy semejantes en cuanto a los cálculos complementarios que deben hacerse. El método de 
descomposición en fracciones simples requiere la búsqueda posterior de cada una de las transformadas 
de las fracciones simples en la tabla aunque estas se puedan reducir a dos con carácter general 
mientras que el método de los residuos da directamente la transformada por aplicación de esta última 
fórmula. Este último método proporciona como contribución a la suma un factor por cada polo sea 
cual sea su orden de multiplicidad, en tanto que los polos de orden "n" necesitan otras tantas 
fracciones simples en el desarrollo de F(p). Ambos métodos deben ser considerados complementarios. 


Anotemos finalmente que puede haber otros tipos de singularidades que no sean polos de 
orden "n" finito como son las singularidades "esenciales" o los puntos "de ramificación”. En ambos 
casos aa establecerse precisiones o modificaciones adicionales en el contorno de Bromwich con 
el fin de extender los resultados anteriores. 


O 2 Las funciones "escalón", "impulso $” y "momento". 
$ 2.1 Función escalón. 
Se puede definir la función escalón unitario u(x) —función "de Heaviside"—en la forma 


1 parax >0 


u() = 


O para x<0 


y su gráfica puede ser representada en la figura 2.1 (a. Esta función u(x) presenta una discontinuidad 
de primera especie en x = O con un salto unidad. La función más general trasladada "a" unidades sobre 
el eje real positivo a > 0 y con una altura "h" se definirá consecuentemente por 


h para x >a 
h u(x - a) = 
0 parax<a 


16 


con la discontinuidad en x = a — figura 2.1 (b — donde hemos multiplicado al escalón unitario par el 
valor real "h". 


y 


figura 2.1 b) 


Esta función será, en principio, enormemente útil para construir y representar 
matemáticamente funciones segmentariamente continuas que son las funciones en las que estamos 
interesados porque son las habituales en prácticamente todos los campos de la técnica. Si queremos 
formalmente construir una función como la representada en la figura 2.2 compuesta por trozos de 
funciones continuas f(x), £(),£() según la definición: 


YG)=0 parax<a, 
YO) = fi) paraa<x<b, Y(x)=f() parab<x<c, Y(x)= fi) parac<x< d, 
Y()=0 parad<x 


figura 2.2 


podemos auxiliarnos de la función escalón definiendo: 
YG) = £16) [ uG—a) - uG—b) ] + £6) [ (16—b) —- uG—O) ] + £360 [ uc) — u(x—d) ] = 
= f1(x) u(x-a) — £1(9 u(«—b)) + £,09 u(—b) — £209 uc) + fa(x) u(x—c) — 309) u(x—d) = 
= £109 uG—a) + [£63)- £69] uG—b) + [f£i69- £,69] uc) — f309) u(x —d). 
estando, en general, este tipo de funciones caracterizado por suma de términos de la forma 
g(x—a) u(x-B) , con g(x) continua en x > 0 
pudiendo ser cero independientemente tanto a. como f. 
La función escalón nos permite representar simplemente un "impulso rectangular” de altura 
"h" y ancho (b—a) en la forma í(x) = h [u(x—a) —u(x—b)]. Está situado en la parte positiva del eje 
centrado en x = = tal y como se representa en la figura 2.3. 
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i(x) = h [u(x-a) - u(x-b)] 


figura 2.3 


Otros impulsos con diferentes formas (triangular, senoidal, etc.) pueden ser formalizados de la 
misma manera. Así podemos, por ejemplo, representar formalmente el tren t(x) de impulsos 
rectangulares de período '"T'" y altura unidad de la figura 2.4 


t(x) = UGT) — u(x-2 T) + u(x-3 T) — u(x-4 T) +... 


que representa un modelo elemental del conjunto múltiple de funciones periódicas cuya representación 
funcional presenta gran interés en la física. 


2T 3T 4T 5T 6T* 


figura 2.4 


Desde el punto de vista físico, la función escalón puede ser vista como una aproximación 
límite de funciones que de forma natural no saltan bruscamente del cero al valor "'h" sino que varían de 
forma continua pero imperceptible a la escala de observación del problema que imponen los demás 
hechos físicos involucrados. Si la variable "x" es la temporal podríamos decir que el salto se verifica 
con una duración menor que la constante de tiempo dominante del sistema. Puede servir de ejemplo 
característico de la función escalón en t = O, el cierre de un circuito con un generador de corriente 
continua. 


El uso de la función escalón nos permitirá disponer de una representación funcional para leyes 
de cargas definidas en ciertos sectores de un intervalo mediantes funciones continuas con la forma 
representada en la figura 2.2. 


8 2.2 Función impulso 6 o función "delta" de Dirac. 


Definiremos ahora la función 
O paraxx*0 
(x) = 
dl 869) dx=1 si 0en. 
que podemos contemplar, si se quiere, como una función simbólica obtenida por paso al límite de una 
función impulso rectangular centrado en x = 0 al ir estrechando su base y aumentando la altura pero 


conservando el valor unidad para el área que intercepta con el eje OX. Esta "génesis" gráfica puede 
visualizarse en la figura 2.5. 
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figura 2.5 


La función más general trasladada "a" unidades sobre el eje real positivo a > 0 y con una altura 
"H" se definirá consecuentemente por 


O para xx 0. 
H S(x-a) = 
JE H5SG)dx=H si 0en. 


La forma de esta función simbólica antes de su idealización es indiferente, puede ser 
consecuencia de un impulso rectangular como en la figura 2.5, de una onda senoidal, de una función 
de Gauss, etc. pero sobreviene al ir estrechando la base y aumentando la altura conservando el valor 
del área "H". El empleo de la función "delta" fue sistematizado por el físico inglés Paul Dirac y 
manifiesta un hecho físico idealizado, enormemente fecundo en el campo de la física y la mecánica. 
No estamos interesados aquí en extendernos acerca del emplazamiento matemático de tales funciones 
que amplían el concepto clásico de función y que tienen un contexto propio en el campo puramente 
matemático. 


En el contexto de la flexión de barras que trataremos más adelante, una carga concentrada es, 
básicamente, la idealización de una carga vertical notable que actúa sobre un espacio muy estrecho 
reducido, en el límite, a un punto y matemáticamente, el objeto que mejor describe este tipo de función 
es justamente el impulso de Dirac o función "9" . Podemos imaginar, por ejemplo una carga uniforme 
Q/e, —altura de un impulso rectangular— actuando sobre un intervalo (a — e/2, a + e/2) que vamos 
reduciendo progresivamente conservando el área que subtiende. Particularmente, en este contexto, lo 
importante será traducir esta génesis en la discontinuidad que en el esfuerzo cortante proporciona una 
carga puntual concentrada y su relación con el resto de las funciones de estado. 


Entre las muchas propiedades que se pueden describir con la función impulso "9" así definida 
destacamos el denominado "teorema de la integral": 


"Para a finito, a > O y siendo f(x) una función continua en x = a y de valor f(a) se verifica 
S E (0) 8 Ga) dx = f(a)" 


Su comprobación es sencilla y representa una importante conclusión para el empleo y manejo de esta 
función. Puesto que el intervalo de definición de la integral de la transformada £ es el [0,00), el 
resultado puede escribirse simplemente f(x) 0(x — a) = f(a) 6(x — a) entendiendo que es en el marco de 
las integraciones donde este resultado debe de ser manejado y sabiendo que $(x — a) con a > O deberá 
ser "aplicado" a funciones continuas en x = a. Es dentro de estas condiciones, en un contexto 
simbólico, en que deben comprenderse resultados tales como: 


9460)=5-x) ;x6(x)=0 ; x S(x-a)=aó(x-a) ; O(ax)= - S(x)  ...etc. 


$8 2.3 Relaciones formales entre las funciones escalón e impulso. Función "momento". 


Según lo dicho anteriormente, el valor de la integral 
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le S(x-a)dx,cona>0 


será nulo para x < a ya que en tal caso el intervalo [0, x] no contendrá a "a" y valdrá la unidad en el 
caso x > a debido a que "a" quedará dentro del intervalo de integración. El resultado de la integral 
coincide con el de la función u(x—a), así que podemos escribir: 


u(x —a) =f, Sx — aJdx 


br d ad 

y admitir que de u(x—a) = S(x— a), donde el operador derivación adquiere un valor puramente formal 
ya que el límite que implica no puede ser definido sobre una discontinuidad. Obsérvese, además, que 
la derivada de funciones tales como las de las gráficas 2.6 (a tienen las formas representadas en la 
figura 2.6 (b. Se trata de aceptar que esta situación se conserva en el límite en el que las funciones se 
convierten en un salto y un impulso ó respectivamente. En la misma línea podemos aceptar este 
A l : ' d 
proceso de derivación formal y considerar la derivada S' del impulso de Dirac em (xa) = S'G— a). 
Las gráficas (c de la figura 2.6 y su paso al límite ocuparían el mismo papel que las derivadas de las 
figuras anteriores. 


f(x) 00 100 


(b 


09) fx) f (9) 


figura 2.6 


Si se quiere proceder a la "construcción" formal de la S' bastará considerar un impulso 
rectangular de ancho e > O centrado en x=a . Esta función responderá a la expresión 


ulx-(a — le ulx-(a+3)] 


si conservamos ahora el área encerrada por el bucle rectangular de manera que sea igual a la unidad 
estrechando la base "s" esta función deberá tender al impulso (x —a). Por lo tanto podremos escribir 


S(x —a) = limo-»o (ulx (a 2) u[x (a+2)D, 


siempre que Jen Cu[x (a 2)] u[x (a + 3) =1 


Considerando derivadas formales a ambos lados de la expresión y conmutando el operador derivada 
con el de paso al límite podremos escribir 


5" (xa) = lime so(5[G: (a- 5)1-ap (a +2) 
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esta nueva "función" puede verse generada por la serie de gráficas de la figura 2.7 en la que se 
representa un doble bucle de "onda" rectangular llevada al límite de dos impulsos 8. Una tal función 6" 


(x—a) además de ser nula para todo x x a, verifica 8 5'(x—a) dx = 0 para cualquier intervalo "n” del 
semieje positivo contenga o no a la abscisa x = a. Sin embargo lo que hace realmente interesante a 
esta función es la particularidad de no ser cero la integral S, x 5 '(x—a) dx cuyo valor es la unidad 


con signo negativo cuando aen, o sea, la de ser 


DE S'(x—aJdx =-— 1 para cualquier e > 0. 


e>0 2>0 


figura 2.7 


Para comprobar este resultado y, al mismo tiempo, conocer mejor la metodología operacional asociada 
a estas funciones trataremos de hallar el lim¿>0 f Ela f() 5'(x-a) dx para cualquier e>0. Efectuando 
E>0 
la integración por partes 
limo [%* fG9 3'(x-a) dx = lim (fo S(x -a)|% +E 5% 60 8(x-a) dx =-f(a) 
go  Y€ SE a—-e a-e 


con f(x) continua en x = a y de valor fa). El primer factor del segundo miembro se anula. Puede 
ahora extenderse el intervalo de esta integral contando con el valor nulo de 0' (x—a) en x z% a, así que: 
Dad f695'(x-—a) dx =-—f(a) , para cualquier e > O, y por tanto 

1 f605 '(x-a) dx =-f (a) ,a>0 

Este sería el correspondiente "teorema de la integral" para la función S'. Ambos "teoremas" pueden 
generalizarse y contemplarse conjuntamente 


0 si azn 

E f60 8 (u-a)dx = 
f(a)si aen , f(x) continua enx=a 
0 si agn 

S, 60 8' (-a)dx = 


—f' (a) si aen , f'(x) continua enx=a 
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que proporcionan el sentido operacional asignado a ambas funciones en el marco de las integraciones. 
Es en este mismo marco es en el que se pueden escribir expresiones tales como 


x6(x) =-05(), 0%) =-090'(-x), ...etc. 
Particularmente se obtiene el resultado ya citado 
Í. x8'(x-a) dx=-1 si aen 


y esto es lo que permite considerar a 9 (x — a) como la "función momento unidad" toda vez que el 
valor absoluto de su momento de primer orden es el valor de la "altura" M que proporcionemos a la 
derivada del impulso 9 


S, Mx58'(x-a) dx =- M si aen 


Por lo dicho, esta "función" nos va a proporcionar una tercera forma usual de solicitación 
sobre los modelos matemáticos de las leyes de cargas en la flexión simple de vigas. Serán los 
momentos de flexión aislados de magnitud "M". 

O 3 Transformadas de las funciones segmentariamente continuas y de los impulsos "5" 
y momentos "3" ". 


8 3.1 Transformadas de la función escalón, impulso e impulso de Dirac. 


La transformada del escalón unitario u(x) será 


>: CONE NN 2 1d px _ 
£ [u(9] =$, eP*Xu (o dx=f, end = 0% Mr 


y la del escalón de altura "q" será £ [q u(x)] = S 


Teniendo en cuenta ahora la traslación en el dominio de la "x" obtenemos 

£ [qu(x -a)] =q€e ”* S 

así que la transformada de un impulso rectangular de altura "h" en (a,b) será 
£ [hi(a,b) ] = £[ h luG—a) — u(x —b)] ] = . (PPP) 

y, consecuentemente, para el definido en (0, b) 


£ [hi(0,b) ]= Ser) 


Para obtener la transformada de 6 bastará considerar su "teorema de la integral" 
£[ Sa) ]= 15 ePX5(x-a)dx=e>”? 
por lo que £ [S(x)] = 1, siendo en general para Qe KR 
¡23009 A dl Q ; Odds cotas Qe *P 


A estos mismos resultados puede llegarse desde la posible génesis formal del impulso ó tanto 
en x = 0 como en x = a considerando el límite un impulso rectangular. 


Análogamente para obtener la transformada de S' (x — a) basta asimismo utilizar su "teorema 


de la integral" y obtenemos directamente 
£ "(x= = E DPXS'(x— dx = — (- ds = map 
[S' («a)] = $, eP*5'(ra) dx=-(pe anses pe 


y, por tanto, £ [0' (x) ] = p. Los resultados generales pueden ser añadidos en la forma 
MO (ii das Mp a MÓO(X—Aeoconccccn... iia Mpe?*” 
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8 3.2 Transformada de f(x —a) u(x-a) y de f(x) u(— a). 
Siendo £ [f (x)] = F(p), en el apartado $ 1.3.2 obtuvimos el resultado 
£ [f(x -—a)]=e P? F(p) , a>0, siendo f(x — a) =0 parax<a 
Es ahora cuando podemos formalizar este mismo resultado estableciendo 
£ [f£ (x- a) uG— a)] = e ?*?P F(p) 


resultado al que también puede llegarse aplicando la definición y realizando un adecuado cambio de 
variable. Cuando a = 0 proporciona el resultado trivial: £ [f (x) uy] = F(p). 


En los cálculos ordinarios, si f(x) es transformable podremos prescindir normalmente de la 
función u(x), es así como las anteriores transformadas obtenidas para las funciones f(x) son válidas 
para las expresiones f(x) u(x). No obstante, el manejo de la función escalón no siempre resulta formal. 
Es el caso de la obtención de transformadas de funciones tales como f(x) u(x — a). Para este tipo de 
transformadas siempre puede acudirse a la definición y utilizar reiteradamente la integración por 
partes pero es mucho más directo utilizar 


£ [£ (9) ubG— a)] = e *PE [f (x + a) uG)] 


que es otra manera de considerar el resultado £ [f (x — a) u(x-a)] =e *P£ [f (x) uo]. 


Las transformadas: 


E A LT DEI e P£ [f£ Gol. 
EEUU. ara raso e P£[£(x+a)] 


son fundamentales para manejar las funciones sectorialmente continuas que se utilizan habitualmente 
en la técnica. Unos sencillos cálculos puede ilustrar las diferentes posibilidades de transformación. 
Sabiendo, por ejemplo, que 
2 1 c 
£11=%,£lx1=35% y £ lc] == 
=> £lx]= y £lel => 


siendo "c'" un valor real conocemos inmediatamente las transformadas: 
E qee. 5 
£ [G-a? u(x -a) = 538%, £[(a) uba)] = 7 €"*” y £ [cu(x-a)] = a e 


Para hallar, sin embargo la transformada de x” u(x—a) podemos utilizar la fórmula correspondiente 
para obtener 


a 2 2a, a 
"(244 =) 


£ Lua -a)]=ePE£[(x+a]=e PE£[x+2ax+ta"]=e o 


o bien podemos descomponer x” en la forma 
x= (xa) +2a(x—a)+al 
y de aquí 


£bé u(x —a)] = £[G—a) u(x —a) + 2 a (x-a) u(x-a) + a“ u(x-a)] = E > an 
8 3.3 Transformadas de las funciones periódicas generadas con una función g(x) en el 


primer período. 


Las transformaciones vistas anteriormente nos permiten plantearnos de forma sencilla la 
transformada de las funciones periódicas f(x) de período "IT" definidas en el semieje positivo —figura 
3.1—. Si es g(x) una función que posee transformada definida en [O, T'], la función periódica f(x) 
generada con esta función será 


6) = 869 u6) + € 3D) UD +2 D) u—2 T) +... 
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O T 27 3T 4T 


figura 3.1 
Podemos prescindir del escalón u(x). Siendo G(p) la transformada de g(x) 
F(p) = G(p) + G(p) e ??+G(p) e *P+G(p) e 77P+...=G(p)(1+0 "+0 P+007P 4...) 
y teniendo en cuenta la suma de los términos de una progresión geométrica ilimitada 


E(p) = G(P) + 


Si se considera el tren de impulsos rectangulares de la figura 3.2 con 


-(5)» 
g(x) = u(x) — u (x — >) y su transformada G(p) = S - —=— 


figura 3.2 


1. EP 1 


su transformada será F(p) = Dú-= Tn = 


ans 
p(1+ e P) 


Si el tren de impulsos es, finalmente, el de la figura 3.3 
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figura 3.3 


e Tp e2Tp 
+ 
p 


con g(x) = u() — 2 u(x —T) + uG—2 T), siendo G(p) = = - 2 , la transformada de la onda 


rectangular no rectificada con período "2 T" será 


el 3YP e PP _ 1 Tp 
a y 
elzP4 e DP Pp 


0 4 Aplicaciones de la transformada. 


Algunas de las aplicaciones tradicionales de la transformada £ han sido hoy totalmente 
superadas por la generalización del uso de ordenadores. La existencia de diferentes programas de 
cálculo numérico y simbólico han dejado fuera de utilidad una gran cantidad de métodos de cálculo 
incluidos el de las transformadas que aquí se han descrito o, por ejemplo, la aplicación de la 
transformada al cálculo de integrales. Resulta sin embargo interesante aprovechar su metodología para 
resolver problemas particulares que requieran un tratamiento sistemático en el campo de las 
ecuaciones diferenciales ordinarias lineales o su ampliación, en algunos casos, a las ecuaciones en 
derivadas parciales y, por extensión, a los sistemas de ecuaciones diferenciales o integro-diferenciales 
simultáneas muy usuales en mecánica, circuitos eléctricos con elementos lineales, problemas 
electromecánicos o sistemas en analogía con los anteriores en cinemática, química etc. 


$ 4.1 Integración de la ecuaciones diferenciales ordinarias lineales con coeficientes 
reales. 


Como se dijo al principio, la transformada de Laplace resulta un potente método de 
simplificación de los cálculos para obtener soluciones particulares de ecuaciones lineales con 
condiciones iniciales. Aunque para algunos tipos de ecuaciones diferenciales lineales con coeficientes 
variables puede aún resultar ventajoso el uso de los métodos de transformación, según el tipo de 
coeficientes, la transformada de Laplace es un método de integración simple y eficaz en el caso de las 
ecuaciones con coeficientes constantes que es el que aquí nos interesa. Podemos describir brevemente 
el procedimiento sobre una ecuación general 


Ca YO + ra YA) + C1 y'(%) + Co yGD) = 1(9) 
Teniendo en cuenta las transformadas de las derivadas y el carácter lineal de £ se llega a 


O 
> (Cr PO Fics +C3p+ C2) y'(0) A Ca y=1(0) 
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si se conocen las condiciones iniciales y(0), y '(0),...... y“"1(0) y siendo £ [f(x)] =F(p), obtendremos 
con simples operaciones aritméticas una función racional, en general irreducible 


£ [y] = v(p) 


CAP" + Cp1p91 4 ...+ 074p +Cp 
y de aquí mediante la transformada inversa 


y = ql [ v(p) 


CpPR + Cp-1 Pr 4+0>...+ 01 P+HC0g 


Describiremos un sencillo ejemplo ilustrativo como es el de hallar la solución y(x) de la 
ecuación diferencial: 


E o x” cosax que satisfaga y(0) = y'(0) = 0 


El procedimiento clásico pasaría por hallar la solución general de la ecuación y ensayar una solución 
particular sugerida por la forma del segundo miembro —método de los coeficientes indeterminados— o 
bien por el método general de variación de las constantes en el que tendríamos que calcular integrales 
del tipo f x? eV3ax cos ax dx y fx? eV3ax cos ax dx. Utilizando los métodos de transformación, 
aplicaremos la transformada a la ecuación 


£ [y"-34y]=£[y")-3a*£[y]=p*£ [y] — p y(0) - y(0)- 3 ad? £ ly] = £ [5 x' cos ax] 
Para hallar la transformada del segundo miembro bastará considerar la transformada de "cos ax” y la 
transformada del producto por x”, obteniendo 


UN 2 3_ 2 
(p?-3a>) £ [y] =20 3a%) así que £ [y] =¿ p a - p 1Pp 3a“p 


(p?+a2)? , p?+a2)? 4 a? (p?+a?)? 4 a? (p?+a?2)? 


y por lo tanto, la solución particular buscada será 


(x) = 2 x senax — — x? cos ax 
y! 8a3 8a2 


En la aplicación de estos métodos al análisis de vigas de rigidez constante con diferentes 
condiciones de sustentación en los extremos, resolveremos la ecuación y* = f(x) para leyes de cargas 
muy generales y este no es un problema de condiciones iniciales sino de condiciones en los extremos o 
"de contorno” por lo que tendremos que habilitar un sistema para adaptar el método al problema 
particular sin perder las notorias ventajas del procedimiento. 


O 5 Análisis de vigas de rigidez constante. 


Consideramos una barra de directriz recta y sección uniforme cargada verticalmente y 
admitimos las hipótesis de existencia de fibras neutras y la conservación de las secciones planas. El 
análisis de la deformación se puede reducir a los desplazamientos verticales que determinan la curva 
deformada "elástica" y a los giros de las secciones que pueden deducirse de esta curva por simple 
derivación. Se considera aquí que las precisiones complementarias y las condiciones que permiten 
reducir este problema a uno unidimensional son ya suficientemente conocidas y nos centraremos en el 
procedimiento para la obtención bajo esas determinadas condiciones de hallar las funciones de estado 
del proceso de flexión utilizando los métodos de trasformación. 


La curva elástica y(x) que define los desplazamientos verticales viene regida por la ecuación: 


3 
y) = 2 1+y (93 
que muestra simplemente la relación proporcional directa de la curvatura de la curva deformada k(x) 
con el momento flector M(x) siendo constante a lo largo de toda la pieza el coeficiente de rigidez de 
la viga a flexión B = E I. Esta ecuación sigue siendo válida cuando la rigidez B es variable. La 
exigencia de considerar pequeños descensos condicionada al comportamiento elástico hacen adecuada 
la ecuación lineal aproximada que resulta de suponer que y '()? = 0. En nuestro caso, con B = cte 
podemos evitar el conocimiento de M(x) derivando dos veces y obteniendo la ecuación 


26 


y Vx )= == 


siendo E la "ley de cargas" verticales, función definida en [O, L] con origen de abscisas en el 
extremo izquierdo de la barra de longitud "L". 


Si w(x) tiene transformada de Laplace W(p), aplicando la transformada a la ecuación anterior 
obtendremos 


e ly "691 =p*£ 1y69] —p? 300) — p?y(0) - p y"(0) - y""(0) = “2 
y de aquí 


1 W(p) , y0) y'(0) ,y"(0) y”'(0) 
A rre ii 


. 1 ¡W , ; A 
Denominando a £* 1 = z(x), las transformadas inversas proporcionarán 


y69 = Ex) +y00) 169 +y00) x +40 2 ¿170 á 


Prescindiendo de la función u(x) ya que a lo largo de toda la longitud "L" de la viga vale 1, las 
sucesivas derivadas de la expresión anterior serán: 


y) =270)+y0 +y"0) +50 y 
y'9=57"9 +y"0) + y"(0) x 
y") =52"09 + y"(0) 
y0=¿28) 
Junto a los desplazamientos y(x), las ecuaciones anteriores nos definen respectivamente el 
resto de las funciones de estado de la flexión simple que son: 


- Los giros de las secciones en el proceso de flexión 6 (x) = arc tg y(x) = y (x) 
- El momento flector M(x) = B y"(x) 
- El esfuerzo de corte o "cortante" T(x)=B y""(x) 


$ 5.1 Las funciones de deformación y esfuerzo en los seis casos estándar. 


Para sistematizar ahora el método requerimos de las condiciones de sustentación de la viga en 
sus extremos que, en principio, son las reflejadas en el siguiente cuadro 2 


SUSTENTACIÓN extremo izquierdo extremo derecho 

DÁ 00 | y '(0)=0 y "(L)=0 y"'1)=0 
> y'(0)=0 y'""(0)=0| YU)=0 | y'(1)=0 

Fi A y(0) =0 y "(0)=0 y(L) = 0 y (L)=0 

Y y(0) =0 y (0) = () y(L) =0 y (L) = () 
S [YO=0 | y'"(0)=0 YU)=0 | y'(L)=0 

7 NS | YO=0 | y'(0)=0 | YU)=0 | y'()=0 

CUADRO 2 


donde se han considerado las combinaciones mecánicamente viables de las tres condiciones tipo: 
"empotramientos perfectos” a los que nos referiremos simplemente como "empotramientos”, "apoyos 
rígidos" o simplemente "apoyos" y "extremos libres". El primero queda caracterizado por descensos y 
giros nulos en el vínculo, el segundo por la anulación de descensos y flectores y el tercero por la 
anulación de flectores y cortantes. 
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Impongamos estas condiciones de vínculo en cada uno de los casos sobre las funciones de estado 
generales representadas anteriormente. 


$ 5.1.1 Viga empotrada-libre. 


La ménsula empotrada en su extremo izquierdo está solicitada por una ley de cargas w(x) representada 
respecto del origen en el extremo izquierdo empotrado. Considerando las condiciones del cuadro 
anterior se obtendrá para y(x) y sus derivadas segunda y tercera 


y69) = 29 20 4 
y"09=352"09 + y"(0) + y"(0) x 

y")=52"9+y"0) 

que deberán satisfacer las condiciones del cuadro para x= L 

y"(L) =¿Z2"(L) + y"(0) + y”(0) L =0 

y"(L)=32"(L) +y"(0)=0 

y de aquí pueden obtenerse los valores y"(0) e y""(0) 

y"(0)=-¿2"1) 

y"(0)=5(Lz"(L) -2"(L)) 

Así que la ecuación de la curva elástica y la de sus derivadas sucesivas serán en este caso 


Lz”” L -z" L zI" L 
0 Dj _ 0 


nr 
(0) 
OE 
6 


By() =2(x) + - 
By'6) =Z'6) +[L z"(L) -z"(L)] x — 0 Xx 
MG)=By"6)=z"6)+Lz"(L)-z"(L)-z"(L) x 
T(9)=By"09=2"09-2"(L) 

$ 5.1.2 Viga libre—-empotrada. 


La deformada de la ménsula en el extremo derecho solicitada por una ley de cargas w(x) y su primera 
derivada serán en este caso 


yG9 = 5 2(%) + y(0) + y'(0) x 
y'(9) = 5209 + y(0) 
Del sistema: y(L) = : Z(L) + y(0)+y(0)L=0 , y(L)= : z(L)+ y'(0) = 0, obtenemos 


y(O=-=¿2b)  ,  y0)=5-xL)+LZ0)] 

por lo que la elástica y sus derivadas sucesivas vendrán expresadas por 

By) =zZ6) + Lz(L) -z(L)-z(L) x 

B y) =7(x) - 2 (L) 

M(9 =B y'"(x) = z"(x) 

T()=By"69=2"(x) 

$ 5.1.3 Viga apoyada—apoyada. 

Para una ley de cargas con origen en su extremo izquierdo w(x) obtendremos para y(x) 
yO) = 5209 +y 0) + LO y 

Del sistema: = XL) + y (0) L + 20 L=0s : z"(L) + y”"(0) L= 0, obtenemos 


o iZ (E iy — 1 pLz”(L L 
yo) =- ¿2 y yO) = 2-22] 
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y la elástica y sus derivadas serán en este caso 


B y(x) = E Ax) q Lz' O - Z a ] a = 2 y 
A Lo pa z(L) 2") 2 
E 


MG) = B y"G) = Na E a x 
T()=B y") =z"Q9 — 


$ 5.1.4 Viga ado 


Se obtiene en este caso 


y(x) Za : Z(x) + Y - x? + y — xi 


y del sistema: E z(L) + o ¡ÉS ¿420 A 1=0, =7"(L)+ y"(0) + y”""(0) L= 0, obtenemos 


Ñ z' 4 3Z(L) nd 1 ma 3z''(L) 
0 =-:2-%2) LOS Tr? 


y con estos jas la deformada y sus le serán 


Bro = 20H EL ED 
By0)= 00) ED Lx +( e =D) y 
z' o PA DA 

L2 L3 2É 
220 3z"(L) 
2L 


MG) = By"(9 ==") +52 
T(9 = By"69 =2"09 + 2 


8 5.1.5 Viga PON 
La curva elástica debe ser ahora 


y) = 529 +y0) +20 y 
A 05 , , A 

L?=0, 2z(L) + y(0) +22 1? =0, obtenemos 
1,3z(L) 3 z'(L) 


'(L)  3z(L) "mu 
yO: ED yO E E) 


y del sistema: + =(L) + y (0) L A 


y con estos valores las funciones de estado sat definidas por 


iS '(L)  3z(1) (1) '(L) 
BI) OA AGA E 
e z'(L)  3z(L) 3z(L)  32'(L) 
BIO) =2lO+ SA a a 
MON — 1 3 e 3z'(L) 
MG) =B y") =2"(9) + ( 2) 
nr pd UA 3 2200 3 "1 
T()=By"6)=7"0) + 42- => 


$ 5.1.6 Viga A 
En este caso la curva elástica debe ser 


y(x) = is = Ux) A y! O y Y == xXx 


Las ecuaciones para hallar y "(0) e y 


(0) serán 
UL) + 20 Té 2 L7=0 , Ez (L)+y"(0)L+ — L? =0, obteniendo 


1 22 E 6z(L) 12 mn 6z'(L) 


mo - EL ED yo EL E, 


y con estos a podemos escribir 
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B y60) =Z(0) + ( = A Ml A 


L2 L3 L2 
MA 3! 2z'(L) 6 z(L) 6z(L) _ 321) 
BO = 10 + (A - Ain E 
Z MY 9 2z'(L) 6 z(L) 12z(L) _ 621) 
MOB a Ap) 


12z(L) — 6z'(L) 


T09 =By"09 =2"9 +2 - E 


Estos resultados nos permiten reducir el estudio de cada uno de los casos a la obtención de la 
función z(x) y a sus tres primeras derivadas y las correspondientes particularizaciones de dos de ellas 
para x = L = longitud de la viga. 


El cuadro 3 resume las funciones de deformación y esfuerzo para cada uno de los casos 
considerados. 


$ 5.2 La función z(x) y sus derivadas. 


Siendo w(x) la ley de cargas, es W(p) = £ [wG9] y z(O) = £” e ] 


En los modelos que habitualmente contemplamos las leyes de cargas están compuestas por sumas de 
tres tipos de acciones en el intervalo [0, L]: 
- Una función segmentariamente continúa en [O, L] con un número finito de saltos 
finitos. 
- Un número finito de cargas Q concentradas. 
- Un número finito de momentos M de flexión aislados. 


cada una de ellas aporta diferentes contribuciones a la función z(x) y a sus derivadas. 


La primera solicitación proporciona a la ley de cargas, en general, suma de términos de la 
forma f(x) u(x — a) o bien f(x —a) u(x — a) con f(x) una función continua siendo la función escalón la 
encargada de formalizar las discontinuidades. Las transformadas de tales funciones se han descrito en 


W(p) 
p* 


el 8 3.2. Básicamente habremos de considerar las transformadas inversas de la función tal y como 


se ha descrito en el 0 5. Debido a la utilización habitual de sectores solicitados por factores de tipo 
polinómico —tanto si son el resultado de aproximaciones de otras funciones más complejas como si 
son acciones directas de segmentos de cargas continuas de envolvente lineal, parabólica o cúbica— 
puede ser muy útil disponer del cuadro 4 dónde se reflejan algunos de los términos usuales de w(x) 
que contribuyen directamente a la función z(x). En dicho cuadro se ha añadido formalmente la función 
u(x) aunque en la mayor parte de casos prácticos prescindiremos de su inclusión —por ser de valor 
unidad a lo largo del intervalo físico en el que trabajamos— dejando solamente como representativo el 
factor que la multiplica. 


Estas aportaciones directas a la función z(x) se puede completar, cuando sea conveniente, con 
algunas sencillas operaciones aritméticas para su más fácil aprovechamiento. Por ejemplo cuando 
tengamos una expresión como la q (x —b) u(x — a) podemos transformarla previamente de la siguiente 
manera: 


q (x —b) u(x — a) = q x u(x — a)- q bu (x-a) = q (x — a) u(x— a) + q (a — b) u(x-— a) 


La segunda solicitación proporciona términos en general de la forma Q S(x —a) y la tercera 


términos M 8(x —a) por lo que los cálculos propios de la transformación se reducen y simplifican 
mucho. A este respecto es interesante destacar las contribuciones a la función z(x) y sus derivadas de 
tres solicitaciones tipo que se representan en el cuadro 5. Como ya hemos visto la derivación de una 
función escalón genera un impulso. Esto se puede poner de manifiesto cuando la viga está solicitada 
por un momento aislado y en ese caso la derivada tercera de z(x) lo contiene adquiriendo solamente un 
sentido formal para posibilitar el que se verifique la igualdad z'“(x) = w(x), no aportando ninguna 
información relevante a los esfuerzos de corte. En el cuadro 5 se ha conservado este formalismo pero 
en los casos prácticos, puede prescindirse de su aportación a la función z""(x). 
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Viga empotrada - libre 


La "L)> z"(L) e E z""(L) e 

2 6 
By) =2'09) + [Lz"(L) -2"(L)] x- 
MG) =B y"(x) =Z"() + Lz"(L) — ZM(L) Z"(L)x 
T()=By"(9) =2"(%) -2"(L) 


B y) = z(x) + 
ea ye 


Viga libre - empotrada 


B y09 = 29 HD - 2D]: 200 


6L 

y 12 E 0) zD). 2 
B yy =2(9 52-22-22, 
M(x) = B y"(x) = e x 
TG) = = B y" (x) = mn aX ) - z! Ea 
Viga apoyada - apoyada 
By69) =209 HD - 207,20 

_ 12 O Z(D) 202 
By) =7100) + NETAS 
MG) =B y'09 =2'69 - =n 
TG) = B y" (x) = => z"(x) — z! 14) 
Viga empotrada - apoyada 

'() _ 3z(L (1) (1) 

By) =209H 2-20 24 2-2 


2 12 


By(=79) + (40 z! En - 2D) + +( 20 30) y 
M0) =8 y'09 2109 + ZO ón a 
T(9 =By"09 =2"09 +1 20 O 

Viga apoyada - TA 

By6)=29 + (2-20), , 20-50), 

By (9) =2109 + q0 - E DA 

MG) =B y") =2"69 + (2D 20, 

T(9 =By"09 =2"09 +22 ao Pen 

Viga empotrada - E 

By) = 29 2 y A 
By0)=209 +2 - ED) + (72 220, 
M0) =8 y'09 =2109 +40 0-40, (20 - En, 
109 =By"09 =2"9 +24 20) - En 


CUADRO 3 


Solicitaciones W(»)p" z(x) 

que) a Go 

q xu(x) ss - Xu) 

qx u() Fl — 109) 

qx u(x) S = Xx uo) 

q u(x a) A a» a ua) 
qxu(x- a) E ea» a ua) 
qx u(x-a) > a a ua) 
qx u(x-a) a ap a) u—a) 

CUADRO 4 
Solicitaciones w(x) y W(p) z(x) y sus derivadas 
q 2(x) => [6-a)' ua) — (——b)' u(—b)] 
ñ———KÁ bea e do 269 = 7 [G Bd a) —(x py 0 b)] 
: 26) =2 (a) uta) - (DJ u(b)] 


segmento de carga 
uniforme zx) = q [G-a) u(x-a) — (——b) u(x—b)] 


dl z(x) = Q En u(x—a) 
e w(x) = Q S(x-a) (a)? 


x=0 a x=L z()=0Q +. UA) 


Wi) =Qe* 
carga "concentrada" z"() = Q (a) u(x-a) 
z'"(x) = Q u(x-a) 
M z(x) = M pEaE u(x-a) 
Ye | 
o a X=L w(x) = M5 (x-a) z'(x) =M (xa) u(x-a) 
W(p=Mpe”?” MN — 
momento "aislado" p) dl z (x) = M u(x-a) 
z'"(x) = M S(x-a) 
CUADRO 5 


$8 5.3 Acciones aisladas en los extremos. Las particularizaciones de la función z(x) para 
los valores de x = L. 


En una viga de longitud "L" el conjunto de las funciones de estado están definidas y tienen 
normalmente sentido para valores "x" en el intervalo (0, L). Esto implica, por ejemplo, que la función 
u(x) tomará siempre el valor unidad y es por lo que normalmente evitaremos utilizar la función u(x) 
cuando multiplica a algún otro factor como veremos en los ejemplos propuestos en el O 6. 


Cuando existen acciones aisladas en el extremo izquierdo origen de abscisas tanto la función 
z(x) como sus derivadas y sus particularizaciones para x=L son recogidas sin ninguna dificultad y su 
contribución correctamente evaluada, a través de ellas, sobre las cuatro funciones de estado. Las 
posibles dificultades sobrevienen en las posibles acciones aisladas en el extremo derecho; esto ocurre 
en los extremos libres o simplemente apoyados en x = L —figura 5.1—, toda vez que su inclusión se 


de 


verificaría a la derecha del extremo derecho que queda, trivialmente, fuera del ámbito del intervalo 


que subtiende la propia pieza. 
_A —G 
x=L 


libre x=L 
apoyado 


figura 5.1 


Existen diferentes maneras de solventar este tipo de problemas pero la más simple es la de 
precisar que en este caso es por la derecha de L por la que hay que evaluar las particularizaciones de 
z(x) y de sus derivadas. Es decir: se debe de considerar siempre z(L) = z(L”), z(L)= z(L”), z'(L) = 
Z(LP, z"(L) = z"(L5). Este criterio puede hacerse extensivo con carácter general pero resulta 
suficiente considerarlo en estas condiciones particulares. Tanto en el desarrollo del apartado $ 5.1 
como en el cuadro 3 se deben de considerar implícitas estas precisiones en las correspondientes 
particularizaciones de z(x) o de sus derivadas y no se ha hecho así para no complicar innecesariamente 
la notación y la exposición del método que es, en cualquier caso, muy simple. Tendremos ocasión de 
aplicar este criterio en alguno de los problemas realizados en el 0 6. 


8 5.4 Precisión acerca de los signos implícitos en el análisis de la flexión simple para las 
cargas y momentos. 


Se habrá advertido que implícitamente hemos venido considerando el origen en el extremo 
izquierdo de la viga y el eje OX como la directriz recta de la misma en el sistema de ejes habitual 
"levógiro": eje OX hacia la derecha y el eje OY hacia arriba. Por otro lado, las cargas las 
consideraremos como positivas en el sentido de solicitación más corriente, esto es: actuando de arriba 
hacia abajo —sentido "gravitatorio"—. Con estos convenios los descensos de cada sección de la pieza 
serán positivos y los giros de las secciones dextrógiros resultan positivos en tanto los levógiros son 
negativos. Esto significa que los resultados de las funciones de deformación no coinciden con los 
geométricos resultantes de la deformación física. Por esa razón se utilizan diferentes convenios para 
que se verifique la coincidencia, bien sea establecer en el origen el eje OY hacia abajo o sentido 
"dextrógiro" de ejes o bien introducir un signo menos en determinadas ecuaciones etc. Estos criterios 
se encuentran indistintamente en los distintos tratados de la flexión según las preferencias de cada uno. 


En cuanto al sentido de las funciones de esfuerzo, estos son completamente independientes de 
las de deformación. Adoptado el sentido gravitatorio positivo para las cargas ya sean concentradas o 
distribuidas, deberemos ahora elegir un signo para los momentos aislados. Tomaremos el sentido 
"levógiro" como positivo para estos momentos lo que resulta congruente con el asignado a las cargas. 
De esta manera, a momentos flectores negativos en un punto "x" les corresponderá en el entorno de 
"x" una incurvación de las fibras superiores —a la fibra neutra— contraídas en tanto las inferiores están 
extendidas. Los resultados numéricos obtenidos para las funciones de esfuerzo coinciden con este 
criterio y muestran los resultados con su signo correspondiente aunque es muy habitual representar 
gráficamente los momentos negativos en sentido contrario. 


Es significativo representar las gráficas de deformación y esfuerzo de cargas y momentos 
aislados indistintamente en una viga cualquiera —por ejemplo la doblemente apoyada de las figura 5.2— 
y representar los signos que resultan de aplicar los criterios establecidos independientemente de cómo 
se desee disponer las representaciones gráficas para ampliar su expresividad o para acomodarla a la 
representación acostumbrada. se puede observar, además, que se ha incluido un impulso "S" en la 
gráfica de cortantes cuando la viga está solicitada por un momento aislado. 
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figura 5.2 
O 6 Ejemplos de aplicación 
Ejemplo 1: Viga empotrada-libre. 


6m 


figura 6.1 


Se considera la ménsula de longitud L = 6 m. con las divisiones y cargas de la figura 6.1. La 
ley de cargas vendrá representada por 


w(x) = (+3) lu) — uG—-3)] + Ex+6) [u(x-3) — u(x-6)] + 2 8' (x-3) — 0' (x-6) + 3 9 (x-6) 
prescindiremos de la función u(x) pues vale la unidad a lo largo de la longitud "L", así mismo 
prescindiremos de la contribuciones de carga y momento aislado en el extremo derecho para incluirlos 


a través de las particularizaciones de z'(L) = z'(L+) y z"(L) = z"((L+). Con este criterio y sencillas 
operaciones aritméticas la ley de cargas puede expresarse en la forma 


w(x) = —x + 3 + 3 u(x-3) + 2 0' (x-3) 
Teniendo en cuenta los cuadros 4 y 5, la función z(x) correspondiente y sus derivadas serán 


2(x) == +30 03)  u3) +(3) (13) 
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z(x) =- e 0 > + > (6-3) u(3) + 2 (x-3) u(x=3) 
z"(x) =-— - Xx + > Xx + > (3) u(x-3) + 2 u(x-3) 


z"()=- > Xx +3 x +3 (x-3) u(x-3) 


: ad 67 65 
Los particularizaciones z'(6) = z'"(6+) = e 1 = e z'"(6) = z"(6+) =9+3=12 nos permiten 
67'*(6+)-2''(6+) _ 79 z''(6+) 

4 6 
cuales podemos reconstruir la elástica y el resto de funciones de estado según el mismo cuadro 


B y(x) =- A 3 u(x3) +(x-3) u(x3) +Lx-2 
By(x)=- = Xx + > Xx + > (31 u(x3) + 2 (4-3) u(x-3) + Z x—-6x 
MG) =- Ñ xo + . Xx + Ñ (3) u(3) + 2 u(x-3) + Ñ -12x 


TO) =-4x2 +3 x +3 (1-3) u(-3) - 12 


obtener, siguiendo el modelo del cuadro 3, los factores = 2, con los 


funciones que pueden representarse junto a la ley de solicitaciones en la figura 6.2 en la que se han 
incluido algunos valores significativos de los cortantes, momentos, giros o  deformadas 
particularizados del conjunto de funciones de estado obtenidas. 


- 12 | cortantes 


flectores 189/2 


B giros de las secciones =B y' 


15813/40 


B descensos de las secciones=B y 


figura 6.2 
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Ejemplo 2: Viga libre-empotrada. 
3t 


a parabólica 

parabólica ' Dl 2t 
"a al eN. Mn a 
NN STEREO 


6 m 
figura 6.3 


Consideramos la ménsula libre en su extremo izquierdo, longitud L = 6 m. con las divisiones y 
cargas de la figura 6.3. Las componentes continuas tienen una envolvente parabólica perteneciendo 


z A Al ll A 
cada uno de los sectores a parábolas de ecuaciones c, = A x y = E (x-6)' respectivamente. La ley 
de cargas será 


wGo)= -— 8()+ ES luG9-u(x3)1+ : (x-6y [u(x-3)-u(x-6)1-2 5(x-3)+2 5'(x-3)+2 S5(x-5) 


prescindiendo de u(x) y de u(x-6) y operando los factores de u(x-3), la función resulta 


w(x)= — 9 (x)+ =x 4 (x-3) u(x-3) — 2 0(x-3) + 2 8'(x-3) + 2 6(x-5) 


así que la función z(x) y su derivada primera serán 


1 6 1 5 1 3 2 1 3 
Xx) =-2 a Se 3) u(x-3) ¿O 3 U6-3)+HX-3) u(x 3)+ 7 (1-5) u(x-5) 
z(x) = Ne. ¿A 3 u(3) - (x-3Y u(x-3)+ 2 (x-3) u(x-3)+ (5) u(x-5) 
Los particularizaciones z(6) = — y z(6)= e nos permiten, según el cuadro 3, obtener 
6 z(6) — (6) = 35, —z (6) = — 2, con lo que 
B y(x) =-3X . x 5 Qs 3)" u(x3) = (e 3 ux-3)+(3) u(x-3)+ + > (x 5) u(x5) + 
217 
10 
By) A - (3) u03)-03) u(3)+2 (43) u(3)+ (5) us) - E 


M(x) = -1+ = x : (3) u(x-3)-2 (x-3) u(x-3)+2 u(x-3)+2 (x-5) u(x-5) 
Tx) = : Xx 2 (x-3) u(x-3)-2 u(x-3)+2 u(x-5) 


funciones de estado que pueden representarse junto a la ley de cargas en la figura 6.4 en la que, 
asimismo, figuran algunos de los valores representativos de las funciones de estado. 
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2 


ley de cargas a , 6 


3 


0 cortantes 


3.25 
1.25 


0_;4 flectores 


¡247/10 


3383/30 B giros de las secciones =B y' 


B descensos de las secciones=B y 


figura 6.4 
Ejemplo 3: Viga doblemente apoyada. 


A 
e | 


5m 
figura 6.5 


Consideramos ahora la viga simplemente apoyada en sus dos extremos de longitud L =5 m. 


con los segmentos y cargas que se muestran en la figura 6.5. La solicitación vendrá definida por la 
función 


Leal 


w(x) = - x [u6) — uG—2)] + 4 [u(x—-2) — u(x-5)] +3 (2) -6(x-2) -2 8(x-3) -2 6(x-4) —0 '(x-5) 


operando y prescindiendo de u(x) y de u(x-5) puede escribirse 
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w(9 =5x —5 (0-2) (1-2) + u(x—2) +3 S(x-2) -5(4-2) 2 S(x-3) 2 S(x-4) 55) 
y la función z(x) y sus dos primeras derivadas serán 
2(x) = (2) 02) + 2 (02) u2) +5 (2) u(2) —5 (2) (12) 
o (x3) u(x3) Ñ (x-4) u(x-4) 
z(x) = = Xx” = (2) u(x-2) + - (2) u(x-2) + - (x-2)” u(x-2) —(x-2) u(x-2) 
(3) u(x-3) — (4) u(u4) 
z"(x) = . x— . (2) u(x-2) + o (2) u(x-2) + 3 (x-2) u(x-2) — u(x-2) — 2 (x-3) u(x-3) 
— 2 (x-4) u(x-4) 


Los valores z(5) = z(5+) = =, z'"(5) = z"(5+) = 31 + 1, nos permiten calcular los factores que afectan 


a este modelo según el cuadro 3 


a 1 1319 
¿2(5+) 25) 


075 


Lo 16 
== +4) =-=— 
30 2 (5 ) 15 


y con ellos, podemos escribir las funciones 
By(9) = 1-2) (2) + L (2) u2) +5 (2) (2) - 5 (2) u(2) 
= (e 3) u(x-3) = (e 47 u(x y+ -- d 
By'G) = = e = A) ul2) + - (2) u(—2) + > (x2)' u(—2) (x-2) u(x-2) 
(3) u(3) - (4) urea) + 2-2 
MG) =2x-2 (2) u(1-2) + 5 (1-2) u(x-2) + 3 (x-2) u(x-2) — u(x-2) — 2 (2-3) u(x-3) 
2 (4) u( 4) - Ex 
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T(x) = . x= , (x-2) u(x-2) + , (x-2) u(x-2) + 3 u(x-2) — 2 u(x-3) — 2 u(x-4) Fa 


que, junto a la ley de cargas, quedan representadas en la figura 6.6 
Ejemplo 4: Viga empotrada-apoyada. 
Se trata de la viga empotrada en su extremo izquierdo y simplemente apoyada en el 


derecho de L = 5 m. de longitud con los sectores y las cargas de la figura 6.7. La ley de 
cargas puede escribirse 


3t 34 


a 


2 mt 


5m 


figura 6.7 
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-3.4 cortantes 


1319/75 AN cr 


B giros de las secciones =B y' 


- 1381/75 


0 B descensos de las secciones=B y 


figura 6.6 


wo) = (=x+3) [u(9 — u(3)] + (5 x—>) lu(e-3) — u(1-5)]+2 S(x-3) +2 8'-3) 2 5(4-5) 
Prescindiendo de u(x) y de u(x-5) y operando con el factor de u(x—3) resulta 
w(x) = -x+3 + : (3) u(x-3)+2 5(x-3)+2 0'(x-3) -2 8'(x-5) 


Las funciones z(x) y sus dos primeras derivadas serán 


1 1 1 1 
Ll) == + AY 0) + (3) u(3) + (3) u(-3) 
20) = ¿4 (3) (3) + (63) (3) + 2 (4-3) u(-3) 
z"(x) =-— - xo + - Xx + Z (3) u(x-3) + 2 (x-3) u(x-3) + 2 u(x-3) 
con los valores z(5) = = y z"(5+) = 26-—2=24 calculamos los factores 
Lo 3 _ 487 1 Do _ 2887 
O 00? 2502) 72926) = 5000 


con lo que podemos escribir las funciones 


By) == A UA) A UE a) a 


o ME 
B y'(x) ia 1 x + 1 xXx + e 0 3) u(x 3) + (Os 3) u(x 3) +2 (x=-3) u(x-3)+ 0 x-— 1004 x 


5 487 2887 
MO) == ++ (E 10)+ 2 (3) u=3)+ 2003) += X 
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1 5 2887 
TD) =-5 x+3x + a (x-3) u(x-3) + 2 u(x-3) — 000 


que son las representadas en la figura 6.8 


3 3 


L, 


0 ley de cargas 


1863/500 
0.726 ¿ 


cortantes 
- 2887/1500 


487/100 


flectores sip 


B giros de las secciones =B y' - 17/10 


0 B descensos de las secciones=B y 
figura 6.8 
Ejemplo 5: Viga apoyada-empotrada. 
2t 2t 


parabólica 


IN 
CU 
N Jl ll M 


figura 6.9 


Se considera la viga de L = 6 m. de longitud definida en [0,6] cargada en su componente 
continua con un sector lineal en [0,2] y otro de envolvente parabólica cuyo máximo de valor 2t. se 


» eE Ll a .. 
encuentra en x = 4 pudiendo ser representada por la ecuación c = — A (x-4)' + 2. Esta solicitación y el 
resto de cargas y momentos aislados vienen representados en la figura 6.9. 


Para esta pieza la ley de cargas vendrá representada por 
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w(x) = x [u(x) — u(x-2)]+( > (x-4)? + 2) [u(x-2) — u(x-6)] -2 8 (x-2) -8'(x-2) + 3 3 (x-4) 
Prescindiendo de u(x) y de u(x-6) y operando con el factor de u(x-2) resulta 


w(x) = x + ( =(x 2) + (x 2) — 2) u(x-2) -2 6 (x-2) -0(x-2) + 3 8 (x-4) 


con lo que la función z(x) y su primera derivada será 


1 5 1 4 1 3 1 2 
z(x) = TOS — 2 u(x 2) + (8 2) u(x-2) ye 2) u(x-2) 3 2) u(x-2) 7 2) 
u(x-2) + a 7 (4) u(x-4) 

z(x) =Lx 57 LY) u(x-2) +2 2 u(x2) =(x 2) u(x2) — (2) u(x=2) 
—Ax-2) u(x-2) + 5 > (ay u(x-4) 
Los valores z(6) = — y z(6)= e nos permiten obtener los factores 
5 as SE 292 
2216) -7 26) ==, 246) - 5216) = ME 
con los que ES escribir las funciones de estado 
By()= (2) (2) + (2) (2) - E (2) u(a-2) — 3 (2) u(x-2) — 5 
(42) u(x-2) +5 > (4) u(x-4) + = x- E x 


By) =>7X a e 2) u(x ya 2)* u(x2) 5 (e 2)" u(x-2) — (x-2)* u(x-2) 
-(x-2) uba) + + = (ay u(x-4) + 252. e _— Xx 
MG) =1 4 7 2) u(x-2) + 2 E, y 2) u(x2) — (4-2) u(x-2) — 2 (x-2) u(x-2) — u(x-2) + 3 (x-4) 
u(x-4) — -58 
T(x) =32 (2) u(x—2) +3 2Y u(x-2)-2 (2-2) u(x-2) -2 u(x-2)+3 us) E 
senoidal 


que son las representadas en la figura 6.10 
Ejemplo 6: Viga doblemente empotrada. 


ALLI | ME 


figura 6.11 
La barra doblemente empotrada de la figura con una longitud L =8 m. está solicitada por una 
sobrecarga de envolvente senoidal, simétrica respecto del eje vertical que pasa por el punto medio de 
la viga x =4m. donde alcanza su máximo de 3 t. Hay además una carga y momentos de flexión 
aislados en los puntos que se muestran en la figura 6.11. La solicitación completa viene determinada 
por la ley de cargas 


w(x) = 3 sen A x 2 8'(x-2) + 4 5 (x-4) +2 5'(x-6) 


3t 


| 
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e 2 parábola 
0 6 
s ley de cargas 2791/1405 
4.224 


0.558 


0 ! 
- 584/405 -1.442 cortantes 1127/1135 


0 


- 2.55 
232/45 flectores 


B giros de las secciones =B y' 


0 B descensos de las secciones=B y 
figura 6.10 


Ñ T T 1 
La transformada de la componente continua wc = 3 sen gXes W.=37 207 Para encontrar la 
8 


1 A ER 
transformada inversa de eS =3 a HoOn podemos proceder previamente a su descomposición en 
8 
fracciones simples, resultando 
T L _ 1536 1 24 1 1536 1 


8ptp24 0) a pp 0D) apto or op? 


por lo que la función z(x) y su primera derivada resultarán ser 


zZ(x) = 2229 sen A x+ - y (2) u(—2) + : (4* u(x4) + (x-6)' u(x-6) 


q13 
zx) = =— OS cos = x+ E Xx == 2 (x-2) u(=-2) + 2 (x-4) u(x4) + 2 (x-6) u(x-6) 
Los valores z(8) = = + 7 pa z'(8) = ma E +24, nos permiten obtener 
192 5 _ A 4 
5218) — 7 ABS Ss En eo a 


que son los factores de este modelo según el cuadro 3. Las funciones de estado en este caso serán 


B y(x 2 Ea a a (x 2Y u(x 2) +3 (e 47 ul 4) + (x 67 u(x 6) 


3 


192 


A )x DES 


42 


T 


6 12 2 1536 
COST2X+—X 
8 T 


3 


By = 5 2 (x-2) u(x-2) + 2 (x-4) u(x-4) + 2 (x-6) u(x-8) 


713 
Ae 
MG) =- A senix+ Ex 2 02) +4 (x-4) u(x 4) +2 06) + (HF +5)-(+T)x 


T(9 =-HÁ cos Ex += +4 (4) 24 


cuyas gráficas junto a la ley de cargas pueden verse en la figura 6.12 


3 


ley de cargas 9.639 


2 


-2 
cortantes 


17.384 17.384 


flectores 


B giros de las secciones =B y' 


B descensos de las secciones=B y 


figura 6,12 


OD 7 Extensión a las vigas "continuas". 


El método planteado se hace extensivo a las vigas de varios tramos separados por apoyos 
rígidos interiores denominadas "vigas continuas” siempre que la rigidez B se conserve constante a lo 
largo de su longitud total. Considerando el extremo izquierdo de toda la pieza como origen de 
abscisas, bastará substituir cada uno de los apoyos interiores por cargas puntuales desconocidas R; en 
cada uno de los valores x=a, x=b, x=c ...en los que están situados los apoyos interiores —figura 7.1—. 
Se procede después a obtener la función y(x) de los desplazamientos con el mismo método empleado 
en función de los valores R¡ y se pasa a obtener estas reacciones desconocidas con las condiciones 
adicionales y(a)=0, y(b)=0, y(c)=0...Se obtendrá así un sistema de Cramer del que se deducen los 
valores de estas reacciones y con ellos se reconstruye la función y(x) y sus derivadas. 
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x=x; R; 
y(x;) =0 
figura 7.1 
$ 7.1 Reducibles a modelos estándar. 


Cuando los extremos de la viga continua siguen los tipos descritos en el cuadro 3 no hay 
ninguna precisión adicional y basta la equivalencia con el modelo de un solo tramo y la resolución del 
sistema de ecuaciones algebraicas correspondiente. Este proceso se describe en los dos ejemplos 
siguientes 


Ejemplo 1 


(a) figura 7.2 (b) 


Para obtener las funciones de estado de la viga continua empotrada-apoyada de dos tramos 
iguales de 8m. de longitud solicitada como se muestra en la figura 7.2 (a), previamente la reduciremos 
a una viga simple de un solo tramo como se muestra en la figura 7.2(b). La ley de cargas será 


w(x) = u69)-u(x-4) + 2 [u-4)ux-12)] + [u(x-16)] + 3 6(x-2) + R S(x-8) + 3 8(x-14) 
Prescindiendo de u(x) y de u(x-16) y operando los primeros factores se escribirá 
w(x) = 1 + u(x-4) -1u(12) + 3 6(x-2) + R S(x-8) + 3 8(x-14) 


con lo que la función z(x) y sus dos primeras derivadas serán 


z(x) = e x + > (2) u(x-2) + e (x-4)* u(x-4) + e (x-8) u(x-8) - (x-12) u(x-12) 
+5 (14) u(14) 

z(x) = - xo + - (2) u(x-2) + - (x-4) u(x-4) + Ñ (x-8) u(x-8) - (x-12Y u(x-12) 
+3 (14) u(14) 

2") = 5x0 +3 (0-2) ue) +5 (14) ua) + R (1-8) u(-8) — 5 (6-12) (12) 

+ 3 (x-14) u(x-14) 


e particularizaciones z(16) = dio + ze R A z'(16) = 240 + 8 R nos permiten obtener 
" 495 sb 805 11 
2z (16) - 5 1(16) = 2 + Rs eE 116)-2 ¿2"(6) == R 


con lo que la función y(x) quedará 
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By(x) = = Xx + > (2) u(x-2) + e (x-4)* u(x-4) + e (x-8)' u(x-8) - (1-12) u(x-12) 


1 495 3 805 11 
+5 (14) u(x-14) + (7 +3R) te, 0 x 
EE de 4 1 ny Al AY 495 3 2 805 11 3_ 
De la condición y(8) = 57 8 +S(6 2) emi 4) tE +7R)8 +( Pr R)8"=0 


se deduce R = — o y con este valor las funciones de estado serán 
1 1 1 989 1 
B y) = Sl Xx + a (2) u(x-2) + Sl (4 u(x-4) e (8 u(8) Sl (x-12)' u(x-12) 
1 3 249 2 3013 3 
+ 3 (14) u(x-14) + 6 2688 
By'G) = Y + > (2) u(x-2) + - (x-4) u(x-4) = (x-8y u(x-8) - (x-12Y u(x-12) 
+3 (14 u—14) + a a 
2 28 896 
MG9 = 32 +3 (4-2) u(-2) + 5 (1-4) u(4) — 72 (0-8) u(-8) —5 (1-12) u(x-12) 
249 6039 » 
+ 3 (x-14) u(x-14) + 2. me * 
TG) = x + 3 u(x-2) + (x-4) u(x-4) — u(x-8) — (x-12) u(x-12) + 3u(14) — — Xx 


funciones que pueden ser representadas como en la figura 7.3 


989/56 


ley de cargas 


8.27 


2515/448 


0; 
i - 4.72 
- 3013/448 - 9.38 
cortantes 
249/28 733/56 


flectores 


B giros de las secciones =B y' 


- 166/7 


B descensos de las secciones=B y 


figura 7.3 
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Ejemplo 2. 


12t 12t 12t 12t 12t 12t 


6m 6m 6m 


(a) figura 7.4 


Consideramos ahora la viga continua de tres tramos iguales de 6 m de longitud cada uno, 
empotrada en sus dos extremos y solicitada por las cargas y sobrecargas que se representan en la 
figura 7.4 (a) y el modelo equivalente de un tramo en (b). La ley de cargas será aquí 
w(x) = 2 x [u() — u(x-6)] + 2 (x-6) [u(x-6) — u(x-12)] + 2 (x-12) [uG-12) — u(x-18)] 

+ 2 0(x-3) + 2 S(x-9) + 2 8(x-15) + Ry 6(x-6) + Ro 6(x-12) + 6(x-6) + 0(x-12) 


Prescindiendo de u(x) y u(x-18) y operando los primeros factores se expresará por 


w(x) = 2 x- 12 u(x-6) — 12 u(x-12) + 2 0(x-3) + 2 5(x-9) + 2 8(x-15) + Ry 6(x-6) + R, 0(x-12) 
+ 9 (x-6) + 8(x-12) 


La función z(x) y su derivada primera serán 


2) = 0 +2 (3) (3) — 7 00) US) + 7 (6) us) + Y (6) u(xS) 

+2 (9) u(x-9) 3 (1-12) u(x-12) + (12) u(-12) + Y (12) u(x-12) +2 (15) u(x—15) 
zx) == (3) u(3) — 2 (x-6) u(x-6) + (1-6) u(x-6) + Y (a 6) u(—6)+ (9) u(x-9) 

- 2 (x-12)" u(x-12) + (x-12) u(x-12) + — (x-12) u(x-12) + (x-15)' u(x-15) 


los factores 
- 1 5119 4 | E _ 6887 10, 7 
18 ? (18) - 108 z(18) = A z(18) 324 2 (18) = ao a e 


nos permiten escribir la deformada 


By() = +3 (3) (3) - 5 (-6)' u(-6) + 5 (6) u(-6) + a (x-6) u(x-6) 
+ z (x-9) u(x-9) > (1-12) u(x-12) + > (x-12) u(x-12) + = (x-12)" u(x-12) 


1 5119 4 2 6887 10 7 
to (x-15y u(x-15) + o +3 Ry +3R2) Xx + rr! es R>) x 
y las soluciones Ry y Ra de las dos ecuaciones y(6)=0 ; y(12)= 0 proporcionan 
_ 1941 . _ 1859 
Lo 0 


con las cuales podemos rescribir la ecuación de la deformada y sus derivadas 


By() = +5 (3) (3) — 5 (6) u(x-6) + 5 (6) u(e6) - E (6) u(e-6) 


100 
+3 (9) u(-9) — 5 (12) u(-12) + 5 (12) u(-12) - Í2 (x-12) u(x-12) 
+3 (15) u(x-15) + e a 
By) =5X'+ (3) (3) - 2 (1-6) u(-6) + (1-6) u(x-6) - (6) ua) 


100 
+ (x-9) u(x-9) — 2 (x-12) u(x-12) + (12) ut—12) - 2 (x-12) u(x-12) + (15) u(x-15) 


100 
877 631 _> 
50 100 
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MG) = 5 X + 2 (43) u(3) — 6 (1-6)? u(x-6) + u(x-6) - H— (4-6) u(-6)+ 2 (1-9) u(x-9) — 6 
(12) u(-12) + u(-12) - (5-12) u(x-12) +2 (4-15) ut-15) + Z-Ly 


T(9 = x?+ 2 u(x-3) — 12 (4-6) u(x-6) - HA u(a-6)+ 2 u(x-9) — 12 (1-12) u(x-12) 


1859 631 
ia u(x-12) + 2 u(x-15) — rd 


que pueden representarse como en la figura 7.5 


12 12 12 


- 1941/50 l 1 1859/50 
ley de cargas 


25.38 24.56 e 


- 3.62 .44 
- 631/50 -13.44 -12.62 3.62 


877/50 19.89 20-82 cortantes 19.18 1073/50 


flectores 


B giros de las secciones =B y' 


B descensos de las secciones=B y 


figura 7.5 


8 7.2 Vigas continuas con voladizos en los extremos. 


La extensión a las vigas de varios tramos nos llevaría a considerar tres nuevos modelos de 
vigas a las que deberán ser reducidas otros tantos modelos con voladizos en uno de sus extremos o en 
ambos. Esto es: podríamos incorporar a los modelos los de vigas libre-apoyada, apoyada-libre y libre- 
libre en la hipótesis necesaria de que se trata de vigas continuas con apoyos rígidos interiores. La 
forma más simple de resolver estos casos es partir de la ecuación inicial ” válida inicialmente para 
cualquier tipo de viga de rigidez constante escrita en la forma 


By) = 2(x) +B y(0) u(o) + By(0)x +5 By'(0) +¿By"0)x' di 
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donde impondremos las condiciones conocidas del vínculo en el extremo izquierdo obligándola 
posteriormente a que satisfagan las condiciones del extremo derecho a través del sistema de 
ecuaciones algebraicas correspondiente. Dos ejercicios nos permitirá mostrar el procedimiento. 


der de e 


Ejemplo 1. 


1t 


ML 


(a) figura 7.6 (b) 


El modelo de viga continua con las dimensiones y cargas mostradas en la figura 7.6 (a) puede 
ser reducido a una viga de un solo tramo representado en la figura 7.6 (b) que no es un modelo 
estandarizado en el cuadro 3. 


La ley de cargas será 
w(x) = S(x) + u(x) — u(x-2) + 2 [u(x-2) — u(x-10)] + Ry 5-2) + 2 5(x-4) + Ro S(x-6) + 2 S(x-8) 
Prescindiendo de u(x) y de u(x-10) y operando los primeros factores resulta 
w(x) = S(x) + 1 + u(x-2) + R, S(x-2) + 2 0(x-4) + Ro 6(x-6) + 2 0(x-8) 
La correspondiente función z(x) será 
z(x) = Ñ Xx + 5 Xx + = (2) u(x-2) + A (2) u(x-2) + Ñ (x-4y u(x-4) + = (x-6)' u(x-6) 
+5 (8) u(x8) 
Teniendo en cuenta las condiciones de extremo izquierdo la función y(x) y sus dos primeras derivadas 


podrán expresarse por 


B y) = Y + ¿+ + = (a 2) u(x2) Sl (2)! u(x-2) + 7 = (4) u(x-4) E Ue 6) u(x-6) 


+ y cds x0) +By(0) x 

e = 5 Fe E al 7 2) u(x-2) +=(x 2" u(x2) + (x-4) u(x-4) pi 2 6) u(x-6) 
+ (x-8) nó=d) E By (0) 

By") =x? + > Xx” +Ry (x-2) u(x-2) + > (2) u(x2) + 2 (x-4) u(x-4) + Ro (x-6) u(x-6) 
+ 2 (x-8) u(x-8) 

Las cuatro condiciones y(2) = 0, y(6) = 0, y(10) = 0, y"(10) = O nos permiten obtener las 


By0)=%,By0)=-H Ri=-E Re=-E 


con las cuales pueden reconstruirse las funciones de estado 


B y) Si - + e a + = 10 u(x-2) + e 6-2 ul) + : (4) u(4) + a (x-6)* u(x-6) 
+= 2-8) 100 42 _ 2 Xx 


By) == + Ha eS +? Er (2) u(—2) + - (2) u(x-2) + (x-4) u(x-4) + e (x-6y u(x-6) 
+ (x-8)' a - > 
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MG)= x+ > Xx + = (2) u(x-2) + > (x2) u(x-2) + 2 (x-4) u(x-4) + - (x=6) u(x-6) 
+ 2 (x-8) u(x-8) 
TC) = 1+x+ = u6—2) + (2) u(x-2) + 2 u(x-4) + - u(x6) + 2 (x-8) u(x-8) 


que vienen representadas en la figura 7.7 


f 63/8 f 45/4 
ley de cargas 


5.125 


-4.875 cortantes 


-6.125 


flectores 


-13/13 B giros de las secciones =B y' -1313 


B descensos de las secciones=B y 


figura 7.7 
Ejemplo 2. 


a 
A A A 


alle 


¡ | UI UR 
HULL rl) (IN 


pr Le p: 


figura 7.8 (b) 
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El modelo de viga continua con las dimensiones y cargas mostradas en la figura 7.8 (a) puede 
ser reducido a una viga de un solo tramo representado en la misma figura (b) que no es un modelo 
estandarizado en el cuadro 3. 


La ley de cargas es 
w(x) = u(q) — u(x-3) + 2 [u(x-3) — u(x-15)] + [u(x-15) — u(x-18)] + 96) — 0'(x-3) + Ry S(x-3) 
— 0(x-9) + R, 5(x-9) — 0'(x-15) + R3 0(x-15) + 0(x-18) 


Prescindiendo de la función u(x) y de la contribución de la carga en el extremo derecho para incluirla a 
través de las particularización de B y""(18) = B y”"(18+) y operando los primeros factores quedará 
w(x) = 1 + u(x-3) — u(x-15) + 9) — 8x3) + Ry S(x-3) — 9(x-9) + R» 0(x-9) — 0'(x-15) 

+ Ra 0(x-15) 

La correspondiente función z(x) es 

z(x) = - Xx + — Xx + e (x-3)* u(x-3) : (3) u(x-3) + = (x-3) u(x-3) > (x-9Y u(x-9) 

+ (9) u(x-9) — Á (15) u(x—15) — 5 (15) u(x-15) + Y (1-15) u(x-15) 


Teniendo en cuenta las condiciones del extremo izquierdo, la función y(x) y sus derivadas serán 


By) = x + 5 Xx + 5 (3) u(x3) o (x3) u(x-3) + = (x3) u(x3) o (x-9)* u(x-9) 

+ (9) ue-9) - (0-15) u(x—15) - 5 (15) u15) + (1-15) u(-15) +B y(0) +B y'(0) x 
By'G) = =x + - Xx + - (3) u(3) — (x-3) u(3) + a (3) u(3) — (x-9) u(-9) 

+2 (9) u(-9) — 2 (15) u(x—15) — (15) u(-15) + Y (15) u(-15) +B y'(0) 

By") =x+ o Xx + o (x-3) u(x-3) — u(x-3) + Ry (4-3) u(x-3) — u(x-9) + Ro (x-9) u(x-9) 

— (15) u(x—15) — u(«-15) + Rs (6-15) u(«-15) 

By") =1+x+ (3) u(x-3) + Ry uG—3) + Ra u(9) — (24-15) u(G-15) + Ra u(x-15) 


El sistema de cinco ecuaciones ecuaciones y"(18) = 0, y”"(18+) + 1 =0, y(3) = 0, y(9) = 0, y145)=0 
nos permite calcular B y(0), B y'(0), Ri, Ra, R3 


171 A O A PT 
EY)? BS RS a O 


y con estos valores podemos reconstruir las funciones y(x) y sus derivadas 


By(x) = ES PEE = (3) u(x3) > (x3) u(x3) = (3) u(x3) > (x-9)* u(x-9) 


Z (9) u(x 9) = (x-15) u(x-15) — : (x-15) u(x-15) — e (x-15Y u(x-15) + — - = x 
By'G) = 5 + - + - (x-3) u(x-3) — (x-3) u(x-3) = (x-3) u(x-3) — (x-9) u(x-9) 

E (9) u(x-9) - 5 (0-15) u(x-15)  (x-15) u(x-15) - E (15) u(-15) - + 
By") =x+ : Xx + > (x3) u(x-3) — u(x-3) = (3) u(x-3) — u(x-9) - (9) u(x-9) 
-5 (15) u(x—15) — u(«-15) - $ (1-15) u(«-15) 
By"G) = 1 +x+ (3) u(x-3) = u(x-3) — u(x-9) — (x-15) u(x-15) Ni u(x-15) 


y representar sus gráficas en la figura 7.9 


O 8 Extensión a las vigas "compuestas". 


Hasta aquí, hemos venido utilizando implícitamente el principio de que en cada abscisa "x" la 
correspondiente sección de una la constituye una unión "rígida" caracterizada por ser unión entre dos 
sectores de la pieza que comparten la misma deformación y el mismo giro a derecha e izquierda de la 
sección lo que implica que la función y(x) sea una función de clase uno. El posible salto discontinuo 
de la derivada segunda quedaría condicionado a la existencia de momentos de flexión aislados medido 
en escala B y el de la derivada tercera condicionado a su vez a la existencia de cargas puntuales en la 
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figura 7.9 


sección correspondiente. Estas condiciones constituirían las llamadas "condiciones de conjugación de 
la unión rígida". Es la hipótesis utilizada en los problemas de flexión simple convencionales por lo que 
habitualmente no se hace referencia explícita a la misma. 


Podemos considerar ahora el acoplamiento entre dos sectores de viga de igual rigidez que 
comparten igual descenso pero dejen libre el giro de ambas secciones y donde los cortantes de la pieza 
no queden afectados. Este sería el caso de una articulación que puede pasar a denominarse una "rótula" 
ideal. Sin entrar en el análisis de los modelos físicos y de los condicionantes elásticos o constructivos 
que permitan hablar de estas uniones entre dos secciones en cualquier punto de una viga de rigidez 
constante, se pueden tratar de forma análoga las vigas que incluyen estas articulaciones interiores y 
que venimos denominando vigas "compuestas". Como se sabe el número máximo de "rótulas" 
incluidas en las vigas está limitado por el grado de hiperestatismo de la estructura. Una rótula más en 
una pieza convertida ya en isostática impide que toda la estructura o sectores determinados de la 
misma sean capaces de respuesta a flexión y por lo tanto dejan de ser "vigas". Este acoplamiento en el 
cual ambos sectores comparten la misma "flecha" pero que consecuentemente el giro de cada sección a 
ambos lados de la articulación es diferente implica que la curva elástica será una función continua 
simplemente. Al añadir que el momento de flexión debe de anularse y que la ley de cortantes no queda 
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afectada por la rótula se completan las "condiciones de conjugación" que definen este otro tipo de 
uniones. 


Limitado el número de estas articulaciones internas por el supuesto de que nos conservemos 
en el dominio de la flexión simple, podemos resolver las cuatro funciones de estado de vigas 
"compuestas" de rigidez constante con la misma metodología partiendo de la fórmula Y para cada uno 
de los sectores en que las rótulas dividen la pieza. La única condición directa que aporta una 
articulación es inicialmente la de anular la segunda derivada, las dos condiciones restantes son la 
igualdad de descensos y cortantes en los extremos de los dos tramos de viga que enlazan. Los 
siguientes ejemplos muestran este procedimiento. 


Ejemplo 1. 
2t 2t 


1 ml 


S 


figura 8.1 (a) 


La rótula de una viga con las dimensiones y cargas de la figura 8.1 separa dos segmentos de la 
pieza que pueden ser considerados inicialmente como dos vigas de 4 m y 2 m con la articulación en 
uno de sus extremos respectivamente, como se muestra en la figura 8.1 (b) y (c). 


Para la pieza (a) la ley de cargas es 


w19 = (3 x+2) [u(9) — ue-2)] + (3 (2-2) +2) [u(s-2) — u(i-4)] =-3 x+2+u(i-2)+ 


y su función z,(x) NA es 
O DN Ap ye 
AO ez > +S (A 2) u(x2) 
Imponiendo las a del empotramiento en el extremo izquierdo, la función B y(x) y sus 
derivadas serán 
By10) =- == x a x 5 (2) u(x-2) +5 =B GRA Xx +2 -B y1"(0) x? 


By) == 450 62) ub=2) + Byr"(0) x+ 5B yo) 
By1"G) = == xXx +x +5 (x2)y u(x-2) + B y1"(0) + B y,"(0) x 
By") = =5 x+2x+ 6-2) u(-2) + B y,”(0) 


Para la viga (b) de la figura 8.1 la ley de cargas es simplemente w»(x) = 1 y su zo(x) = = x 


con lo que teniendo en cuenta la única condición de su extremo izquierdo y"(0) = 0, la función B y(x) 
y sus derivadas serán 


Bya(x) => x'+B ya(0) + B y2((0) x +5B y2"(0) x? 
By2(9=7 + B y2(0) +5 B y2"(0) x? 

B y") = > x+By"(0) x 

B y2"(x) =x + B y2"1(0) 
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Las cinco condiciones y,"(4) = 0, y2(2) = 0, y2(2) = 0, y1(4) = y2(0), y1"(4) = y2"(0) nos 
permiten obtener 


By"(0)=%, By"(0)=- 32, By0)=32, By(0)=- E, By."(0)=% 
Con estos valores podemos e las funciones de estado en ambos sectores 

B yu(x) =-= Xx ME > al + (-2)' u(—2) EAS 

B y1(x) =-2 Xx +5 + ay u(x—2) + xx 

B y, "(x) = i E +5 (2) u(-2) ACME 


By") = -; x + 2 x + (x-2) u(x-2) — — 


302 73 53 
+ 3 


By 100 
Pa 
BO tes 
By2"0) => x x+a9% 

By") = ao 


Bastará ahora considerar la función B y:(x) con sus derivada para valores de xe[0,4) y B y2(x-4) con 
sus derivadas en x e (4, 6] o, más formalmente, reconstruir las cuatro funciones de estado de toda la 
viga partiendo de la forma 

y(x) = y109 [u(x) — u(4)] + y2(x-4) [u(x-4) — u(x-6)] = y1(x) + (y20—-4) — y 19) u4) 


en la que en los últimos factores hemos prescindido ya de u(x) y de u(x-6). De esta manera las 
deformaciones y esfuerzos de toda la barra quedarán determinados por las funciones 


dl 1 127 302 
BY ==>0 ii (y A a +12 (4) + - L (4) 
53 3, 1 5 1 64 2,127 

ra0 0 ca» nx 2 a o o 73 53 1 1 
A o E 
E 6 2] 3% ol Ju(-4) 

a A A 2 
Ri (2) u(x 2+2 aL ES A xX ¿ 2) E 
E 


dl 127 53 1 127 
Tx) =-5 e +2x + (2-2) u(x-2) - 3 +[x+3 +5 x= 2x- (x-2) +37] u(%4) 


que pueden ser representadas en la figura 8.2 


Ejemplo 2. 


np 


3 t/ml 


ADAN y | 


La viga doblemente empotrada de dos tramos tiene dos articulaciones a 2m de sus extremos, 
tal y como se muestra en la figura 8.3 junto al resto de dimensiones y solicitaciones. Puede 
descomponerse en tres vigas como se muestra en la misma figura. 


figura 8.3 


j áí : 
Para la viga empotrada-articulada la ley de cargas es wi(x) = 1 y su z1(x) = Pm x”. Teniendo en cuenta 
las condiciones de su extremo izquierdo la función y:(x) y sus derivadas serán 
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B y10) = e xi + > B y,"(0) x? + - B y,""(0) x* 

By) = - x+By,"(0) x + > B y,""(0) x? 

By"0)=5 x2+B y,"(0) + B y,"(0) x 

By,"()=x+B y1"(0) 

Para la viga doblemente articulada en sus extremos la función za(x) es 

20) = 3 +2 (3) u(3) + (1-6) u(«6) 

así que la función y2(x) y sus derivadas serán 

By) =3 +2 (3) u(3) + E (6) u(«-6) + B ya(0) + B y2((0) x +2B y2""(0) x? 
By209=7 +2 (3) u03) +5 (6) u(x-6) + B y2(0) +5 B y2"(0) x? 
B y20) =5 2 +3 (3) u(-3) + R (1-6) u(x-6) + B y2""(0) x 

B y") = 3 x + 3 u(-3) + R u(x-6) + B y2”"(0) 


8 
1 
2 
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La viga articulada empotrada la ley de cargas es wa(x) = 1 y su za(x) = = x”. Teniendo en cuenta la 
condición de su extremo izquierdo la función yx(x) y sus derivadas serán 

B ya) = - x + B y3(0) + B y3((0) x + - B y3""(0) 

By309=3 0 +B ya(0) +5 B y3"(0) x? 

By") =% +B ya"(0) x 

B ya") =x + B y3"(0) 

Las ecuaciones 

y1(2) = 0, y2"(12) = 0, y3(2) = 0, y3!(2) = 0 


que son las condiciones del extremo derecho de cada una de las piezas, la condición adicional 
y2(6) = 0 correspondiente al apoyo interior, más las condiciones propias de la articulación 


y1(2) = y(0), y1"(2) = y2""(0), y2(12) = ya(0), y2""(12) = y3""(0) 
son las nueve ecuaciones que nos permiten obtener 


"ran 7781 men 8677 2519 oy _ 17615 mena 6885 2071 

By: (0) = Tze> B y. (Me PO) Or a? PY (Mi Br ==2 
ras 18415 mens 5541 7 11259 

BIA0)E= 3 BE) RSS 


con lo que podemos reconstruir las funciones de estado de cada sector de la viga original considerado 
como una viga independiente 


21 4,17781 2 18677 3 
By) +12 


By.) =: eS 

Byr9=3 +52 x 

By1"(x) = x— E 

Bya)=1 +1 (3 ula) 12% (6) uo) + EP DOE 10 
A 
By) =7 +3 (3) u(3) - 2 (6) uo) + E y 

B y2"(9) = 3 x + 3 u(x-3) — qe59 u(x-6) + 6885 


448 896 
2071 18415 E 15541 3 


Xx 
112 1344 6 896 


1 
By) =h x+ 


1 3 18415 15541 > 
Bys(x)=7 x"- S on 
ya (9) 6 1344 2 89% 
1 > 5541 
Byv3"(x)=7% x*+=—=x 
ya"(x) 2 896 
5541 
B y3"(x) = x + == 
ya"(x) 896 


Pueden ahora considerarse los correspondientes valores relativos de 'x" para obtener cualquiera de la 
funciones de estado en cualquiera de los tres tramos en que se ha descompuesto la viga original o, si se 
quiere, restableciendo el origen en el extremo izquierdo, obtener los valores del primero para abscisas 
xel0,2), los de y2(x-2) y sus derivadas para valores xe(2,14) y los del tercer tramo en la forma 
y3(x-14) y sus derivadas para xe(14,16]. En este caso para x=2 y x = 14 se verifica un salto finito en 
la derivada primera entre los dos primeros tramos y entre el segundo y tercero, coincidiendo el resto 
de los valores de las y(x), y"(x) e y”"(x) por lo que no queda en ningún caso valores indefinidos. 


Finalmente se restablece la unidad de la pieza reconstruyendo la función 
B yQ) = B ya(9 + B [y2a2) — y109] u—2) + [ya(x-14) — ya(x-2)] u(x-14) 
y sus derivadas. 


De cualquier manera pueden realizarse las gráficas de las funciones de estado que se muestran en la 
figura 8.4. 
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0D 9 Las líneas de influencia. 


Nos planteamos ahora la posibilidad de considerar como variable la posición de una o varias 
cargas y como parámetro cualquier abscisa "x" de una viga de longitud "L", Hay un variado conjunto 
de problemas que explotan esta posibilidad. Para simplificar la exposición consideremos el caso de 
una carga puntual concentrada Q actuando sobre una viga empotrada-apoyada con rigidez a la flexión 
constante a lo largo de toda su longitud "L" —figura 9.1— y centrémonos, por ejemplo, en el momento 
flector de la pieza B y"(x) en cualquier abscisa xe[0, L]. Este momento vendrá expresado según el 
cuadro 3 por la fórmula 
z''(L) _ 3z(L) A 3z(L) 3 z''(L) 


MG) =B y") =2"9 442-200, (0-0M,) 
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Para una posición fija y bien determinada x = a de la carga, la particularizaciones de la función z(x) y 
de su derivada segunda serán simplemente z(L) = - Q (L-a)”, z"(L) = Q (L-a) con lo que el momento 
flector será 


M(9) = Q (a) u(x-a) + 5 Q (L-a)- 350 (L-3)'+( 5 Q(L-a) - FQ(L-a))x 


Observamos que el segundo miembro es explícitamente una función de dos variables: la posición de la 
abscisa "x" donde se desea conocer el esfuerzo de flexión y la posición "a" de la carga Q, así que sería 
más propio escribir M(x, a) y no se hace así porque habitualmente se presume que la posición de la 
carga es fija. 


1 
2 L2 


1 
2 L3 


Si invertimos la situación y consideramos la x = x; fija y bien determinada, se comprende que se desea 
obtener este esfuerzo en dicha sección así que dispondremos nuevamente de una función con una 
variable pero esta es ahora la posición de la carga "a" 


M(a) = Q (xa) u(x-a) + + Q (L-a) - 50 (L-a) +( 5 Q(L-a) - 7 Q(L-a))x 


2L2 
Esta es en esencia el concepto de "línea de influencia", en este caso específicamente la línea de 
influencia del momento flector en la sección "x;" al desplazarse una carga concentrada Q que podemos 
notar Mx, (a) = M (x, a) 


1 
2 L3 


En general, la variable "a" puede tomar valores continuos —o discretos— en un subintervalo cualquiera 
de [0, L] o bien a lo largo de toda la viga. En la línea de influencia debe precisarse en primer lugar la 
función de deformación o esfuerzo al que está referida, después la sección en la que está ubicada, la 
carga O la ley de cargas que varían y finalmente el intervalo de variación de estas. 


Esta nueva situación nos obligará a considerar los valores que toma la función escalón u(x-a) 
con "a" variable y "x" fijo, o sea: los valores que toma u(x;-a). El siguiente cuadro 6 precisa dichos 
nz?" 


valores comparándolos con los acostumbrados ( "a" fijo, ''x" variable ) de manera que cuando "a" es 
variable y "x" fijo, u(x—a) toma el valor cero para a > x y el valor uno para a < x. 


"x" variable, "a" fijo "x" fijo, "a" variable 

u(x - a) u(x - a) 

A RA a 4 
AS 2 AA xk A A A 

o a x (0) x a 

u(x - a) u(x - a) 
A RA a re | 
o x a o a X 

CUADRO 6 


Independientemente de las posiciones posibles de las cargas variables en el caso, por ejemplo, 
de empotramientos o apoyo rígidos las líneas de influencia incluyen también a los extremos en 
cualquiera de los casos. De hecho es muy habitual que en la flexión de vigas estemos muy interesados 
en la variación de algunos de los esfuerzos o deformaciones en los vínculos extremos (reacciones o 
momentos en los empotramientos, reacciones en los apoyos, giros en los apoyos...etc.). Si en el 
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ejercicio propuesto hacemos x; = O y a e(0, L] obtenemos la línea de influencia del momento de 
empotramiento cuando la carga "Q" se desplaza a lo largo de toda la viga desde el extremo izquierdo 
al derecho 


Mg(a) = Q+ +Q (L-a)- 7 5Q (La =7 04 -Q2+Qa+Q 
De manera análoga PUR el cuadro 3, o simplemente derivando, el cortante en la pieza será 
T(9 = Q ua) + 35 Q(L-a?- 7 Q(L-a) 


asi que las líneas de ER de las reacciones en los extremos, por ejemplo, serían 


a 0 


149=00 5 QU - QU y Qs de 
— Tia) =- 35 Q(L-a+ 2Q(L-a)= 35 Qa- Qa+Q 


El problema que estamos tratando, sencillo pero relevante, puede ser generalizado para otra 
distribución de cargas más compleja siempre que, finalmente, su variación pueda ser expresada en 
función de una sola variable tal como la "a" del ejemplo o bien, alternando cargas móviles con cargas 
que permanecen fijas que es un ob que se presenta de forma habitual. Estos procedimientos 
permiten una gran cantidad de variaciones con múltiples aplicaciones. 


Se comprende que no solo pueden obtenerse líneas de influencia de esfuerzos que es lo usual 
sino también de las propias deformaciones. En el problema que nos sirve de guía podemos, por 
ejemplo, interesarnos por la variación del giro de la sección en el apoyo extremo cuando la carga 
varía. En este caso, siguiendo el cuadro 3, los giros de cualquier sección vendrán determinados por 


3zx(L) 32z'"(L) 


2 - 3z(L) 
By) = 206 ua (EL ED (O 
así que 
my Y "()  3z(L 3z(L)  32z"(L 
By(L)= 201 O E) (20 20) pp 


con lo que teniendo en cuenta los valores z(L) y z'(L) obtendremos finalmente 
1 3. 3 2 1 2 
= — = - + = 25 
B 0.(a) E Qa ¿Qa QLa ¿QL 
Es fácil, además, extraer ahora de las líneas de influencia diferentes datos de interés general; así, se 
puede deducir de la línea de influencia del momento en el empotramiento My(a) hallada anteriormente 
la posición de la carga "Q” para que este momento sea máximo sin más que obtener las soluciones de 


35 Mo(a)= 0 -FQa+Q=0 


és 3 , Da 
Los resultados de esta ecuación son ( 1 + E ) L, de los cuales pertenece al intervalo [0, L] la solución 


(1- E ) L = 0.422 L que es la posición de la carga que proporciona ese momento máximo. 


Resumiendo: Las líneas de influencia de una distribución de cargas en una determinada 
sección de una viga se definen de una manera general como las variaciones del desplazamiento, giro, 
momento flector o esfuerzo cortante en esta sección cuando la distribución de cargas varía —en función 
de una sola variable— en un sector del intervalo que subtiende la pieza o eventualmente a lo largo de 
toda su longitud. Cuando la sección se fija en los extremos es cuando obtendremos las variaciones de 
las funciones de estado en esos vínculos. Las líneas de influencia simplifican extraordinariamente los 
cálculos cuando se desea obtener resultados sección-posición de las cargas y son importantes en la 
determinación de esfuerzos máximos en determinadas secciones sometidas a test. Es, como se puede 
comprender, un instrumento imprescindible para el análisis de vigas solicitadas por cargas móviles o 
que pueden variar circunstancialmente de posición. Habitualmente se encuentran en los análisis los 
casos de acciones unitarias tanto momentos como cargas o impulsos finitos de sobrecargas. A los 
valores obtenidos para los esfuerzos con acciones unitarias se le han denominado a veces "valores de 
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influencia". Es fácil comprender que al valor de los efectos debido a acciones no unitarias pueden 
deducirse fácilmente de estos valores. 


Dos ejemplos sencillos completarán la descripción de estos cálculos 


Ejemplo 1. 
E EEE CR 


1t 4t 
2 mt 2 mt 
PEE - APP ES: - PE 4 le 4 
HKEE E áÁá á á y á gTá gá— 


(AíA< 2 K2 << <ÁKÁáÁ 
áíKÁE E <<  á áá]á]á á]áH. 


8m 
figura 9.2 8m 


En este ejemplo consideramos la viga de dos tramos empotrada en sus extremos con las 
dimensiones y solicitaciones de la figura 9.2 en la que una carga puntual unitaria está a un distancia 
"a" del empotramiento izquierdo. Se trata de obtener las líneas de influencia de los momentos de 
empotramiento y las reacciones en los extremos cuando la carga unitaria varía en el primer tramo 
ae[0,4). Como aplicaciones complementarias la de encontrar —si las hay— la posición de esta carga 
para que se anule el momento en el empotramiento izquierdo y la posición para que el momento en el 
empotramiento derecho sea máximo. 


La ley de cargas es aquí w(x) = S(x—a) + 2 8(x-4) + R 5(x-4) y la función z(x) y su primera 
derivada serán 


z(x) = - (a) u(x-a) + (x-4) u(x-4) + - R (4) u(x4) 
2(x) = + (ay u(xa) + 2 (1-4) u(-4) + +R (:-4) u(a) 


Los factores 


1, 3 21d 3 1.2, 1 1 1. 1 de 103, 1 2 1 11 
- == = —a a+ tar 2R+"; — =— =- — 4 + —a -— =R-= 
328) ¿48)= 338 1 +79+R+y; 748) ,208) 1530 * 128% 7 7 R-= 39 
son los que entran a formar parte de la función B y(x) que es 
1 1 1 1 1 1 1 
B y(x) = E (a) ua) + (x-4) u(x-4) + ¿R (47 u(4) + ( al a — > a+ ¿AF GRE eS 
1 3 1 2 1 11 3 
a des A > 
Serra ero 


ada P ál 3 en 
La condición y(4) = 0 nos permite obtener R = En ar — 56 a” y con este valor podemos escribir las 
funciones de estado 


B y(x) = - (a) ua) + (x-4)' u(x-4) + - ( = a” — = a) (x-4) u(x-4) 


1 3 1 > 1 1,1 3 3.) 1 2 1 3 1 2 1 11 3 
Gia? a? e Er HA 1596 * tii 
By) = 7 ay u(x-a) + 2 (x-4) u(x-4) ás Js” u(x-4) 

1 3 1 > 1 1,1 3 3 3 1 1 3) 1 2 1 11 2 
+2 ira 162 A e eS 
M() = (<a) u(x«-a) + 2 u(x-4) (E 4) (4) u(x-4) 

1 3 1 > 1 1,1 3 3 3 1 1 3 1 2 1 11 
+2 (78-38 a TE JEFE 1560 + a? E 45) 

==: => 32 2 e a, 3 53 2 a 
T(x) = uba) + (5 A Ts Ju 4 +6l 1536 E 128 ek 2 
Podemos evaluar ahora 

a 3.3 7.2 1 
M0) = Mole) = y a 15? a , con ae[0,4) 

mes IN o 2 2 
M(8) = Ma(a) = La con ae[0,4) 


T(0) = Tola) == a + Zara E con ae[0,4) 
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— T(8) = - Tala) = — == a+ = a+ = con ae[0,4) 


que son las líneas de influencia de los momentos y reacciones respectivamente en los 
empotramientos. 

1 
La función My(a) = iS = = 2d e = 0 admite efectivamente una raíz real en el intervalo de 
variación de "a", su elo numérico es a = 1.141. En las proximidades de ese punto el momento en el 
empotramiento izquierdo es nulo. 


E 1 s e A E E 7 
La ecuación ¿a ¿a= O donde el ae miembro es la derivada de la función Ma(a) tiene una raíz 


"a " 


en el intervalo de variación de que es =a = 6.222 donde el momento tiene el mayor valor absoluto 


que es - 22 - 0.648. 


Ejemplo 2. 
PC PC 
tt 1t 
pi 4 6 4 6 + 
E 10m y figura 9.3 Ñ 10m Ñ 


Aquí tratamos el caso de la viga doblemente empotrada de 10m de longitud con una 
articulación a 4m del extremo izquierdo solicitada por una carga puntual unitaria que varía en el 
intervalo (4,10] entre la rótula y su extremo derecho —figura 9.3—. Se trata de obtener la variación del 
descenso de la articulación en el proceso de variación de la carga ( línea de influencia del 
desplazamiento de la rótula ) y las líneas de influencia de los momentos de empotramiento. Como 
aplicación complementaria nos proponemos hallar la posición de la carga —si la hay— para que ambos 
momentos de empotramiento sean iguales y la ley de flectores correspondiente a ese caso. 
Acompañaremos aquí el cálculo de estas funciones con las correspondientes gráficas. 


Las funciones de estado de los dos sectores segregados por la rótula serán 
By0)=5By,"(0)7+2B y,"(0) x” 
By1(x)=B y1"(0) x + > B y,""(0) x” 
B y1"G) = B y1"(0) + B y1"(0) x 
B y,"(x) = B y1”"(0) 
B ya(x) = 5 (xa)? u(a-a) + B ya(0) + B y21(0) x +2B y2""(0) x” 
B y2(x) = > (xa) ua) + B y2'(0) + > B y,""(0) x” 
B y2"G) = (x-a) u(x-a) + B y2""(0) x 
B y2"(x) = ua) + B y2"(0) 


De las condiciones B y,"(4) = 0, B y2(6) = 0, B y2(6) = 0, B y¡(4) = B y2(0), B y,""(4) = B y2""(0), se 
obtienen fácilmente 


27 108 1 3,027 27 

” A Eo q a mu A $ AS 

B y1"(0) pa 359 hs (0) 560 “14097 35 
144 76 3 1 2 177 144 1 3 27 27 
B y.(0 ria Bw(0M=2a-l2+ 2-2 Bywr0m==a-%a-2 
yA(0) = 35 2 + 3> y21(0) 280 2 70 35 9 Y2 (0) 560 7 1x0 35 


con lo que pueden reconstruirse las funciones de estado anteriores o bien formalizarse el conjunto de 
ambos tramos con las funciones B y(x) = B y/0) + (B y2() — B y1(%) ) u(x-4) y sus derivadas, donde 
se han de terner en cuenta los valores hallados anteriormente. 


Los descensos de la articulación vendrán determinados por la función B y,(4) = B y2(0), o sea 
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Ñ 2 4 3 144 576 
EU) =BN Fi 


que es la línea de influencia del descenso de la articulación y su gráfica será la representada en la 
figura 9.4 


B y (2)=B y,(0) 


16 
1” 
a 
10 
z 
6 
4 
2 
0 =>a 
0 1 2 3 4 5 6 
figura 9.4 
Las líneas de influencia de los dos momentos de empotramiento serán 
. 1 3 27 108 " 3 3,11 48 
= == =—a+ —: = == ¿== + — —Ñ 
B y1"(0) =M;, (8) = 7,7 4358 + 35 5By2(6)= M2,= 395 8 +79 2+ 35 
y sus gráficas se representan en la figura 9.5 
M;,(a) M, (a) 
0 6 
3 L6 
L4 
25 
ua 
2 1 
15 0% 
06 
1 
04 
sé 02 
o >a 0 >a 
0 1 2 3 4 5 6 0 1 2 3 4 5 6 
a 
figura 9.5 


La igualdad de ambos proporcionan dentro del intervalo de variación de la carga además de la 


“e. E E 8 Ez 
solución a = 6, la raíz a = 2 en la que ambos momentos son iguales a ¿ — 1.6. Para esta posición de la 


carga la gráfica de momentos flectores es la que se presenta en la figura 9.6 


figura 9.6 5 
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